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ABSTRACT 

This paper commences a series of investigations into the stellar populations of local elliptical galaxies as de- 
termined from their integrated spectra. The goal of the series is to determine the star formation and chemical 
evolution histories of present-day elliptical galaxies. The primary galaxy sample analyzed is that of Gonzalez 
(1993, G93), which consists of 39 ellipticals drawn primarily from the local field and nearby groups, plus the 
bulge of Messier 3L Single-stellar-population (SSP) equivalent ages, metallicities, and abundance ratios are de- 
rived from H/3, Mgb, and (Fe) line strengths using an extension of the Worthey (1994) models that incorporates 
non-solar line-strength "response functions" by Tripicco & Bell (1995). These functions account for changes in 
the Lick/IDS indices caused by non-solar abundance ratios, allowing us to correct the Worthey (1994) models for 
the enhancements of Mg and other a-like elements relative to the Fe-peak elements. 

SSP-equivalent ages of the G93 ellipticals are found to vary widely, 1.5 ^ f ^ 18 Gyr, while metallicities [Z/H] 
and enhancement ratios, [E/Fe] are strongly peaked around ([Z/H]) = +0.26 and ([E/Fe]) = +0.20 (in an aperture 
of radius r(,/8). The enhancement ratios [E/Fe] are milder than previous estimates, owing to the application of 
non-solar abundance corrections to both Mgb and (Fe) for the first time. While [E/Fe] is usually > 0, it is not the 
"E" elements that are actually enhanced but rather the Fe-peak elements that are depressed; this serves not only 
to weaken (Fe) but also to strengthen Mgb, acco unting for the ove rall generally mild enhancements. Based on 
index strengths from the Lick/IDS galaxy library (Irager et al. 1998), C is not depressed with Fe but rather seems 
to be on a par with other elements such as Mg in the "E" group. Gradients in stellar populations within galaxies 
are found to be mild, with SSP-equivalent age decreasing by 25%, metallicity decreasing by ([Z/H]) = 0.20 dex, 
and [E/Fe] remaining nearly constant out to an aperture of radius r^/l for nearly all systems. 

Our ages have an overall zeropoint uncertainty of at least ^25% due to uncertainties in the stellar evolution 
prescription, the oxygen abundance, the effect of [E/Fe] ^''0 on the isochrones, and other unknowns. However, 
the relative age rankings of stellar populations should be largely unaffected by these errors. In particular, the 
large spread in ages appears to be real and cannot be explained by contamination of H/? by blue stragglers or 
hot horizontal branch stars, or by fill-in of HP by emission. Correlations between these derived SSP-equivalent 
parameters and other galaxy observables will be discussed in future papers. 

Subject headings: galaxies: elliptical and lenticular, cD — galaxies: stellar content — galaxies: abundances — 
galaxies: evolution 



1. EMTRODUCTION 

This paper is the first in a series on the stellar populations 
of local field and group elliptical galaxies based on the high- 
quaUty spectral data of Gonzalez (1993; G93). The present 
paper concentrates on deriving improved stellar population pa- 
rameters by correcting existing population models for the ef- 
fects of non-solar abundance ratios. The major roadblock to 
population synthesis models of elliptical galaxies is the fact 
that the effects of age and metallicity are nearly degenerate 



in the spectra o f old stellar populations ( Faber 1972 1973; 



O'Connell 1980; Rose 1985 



Renzini 1986). However, it was 



early noted that certain spectral features are more sensitive to 
age than metallicity (e.g., the Ba lmer lines [ O'Connell 1980| ; 
Rabin 19821 ; [Burstein et al. 1984| ; |Rose"l985p , and Sr II A4077 
[ [Rose 19'^ ]), and hope grew that such featiires might be able 
to break the degeneracy if accurately calibrated. (At about the 
same time, several workers were also using Balmer lines to 
discover strong bursts of star formation in so-called "E+A" or 



"post-starburst" galaxies [pressler and Gunn 1983 ; Couch & 



'Carnegie Starr Fellow 
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Sharpies 1987; [Schweizer et al. 199Cl| ], but these applications 
always implicitly assumed solar metallicity.) 

Our ability to decouple age and metallicity in integrated 
spectra has greatly improved over the last decade, due to 
three developments. In the late 1980's, interior models of 
super-solar-metalli c ity stellar evolution becam e available (e.g., 
VandenBerg 1985L [VandenBerg & Bell 19851 V andenBerg & 
Laskarides 1987; pertelU et al. 1994D . Next, the Lick/IDS stel- 
lar absorption-line survey provided empirical polynomial fitting 
functions for a set of standardized absorption-line indices as a 
function of st ellar temperature, gravity, and metallicity (G or- 
gas et al. 1993; |Worthey et al. 1994| ). Finally, an extensive grid 
of theoretical model atmospheres and stellar flux distributions 
was provided by Kurucz (1992) for stars over a wide range of 
temperatures and metallicities. With these three ingredients, 
it finally became possible to compute absorption-line strengths 
from first principles for single -burst stellar po pulations (SSPs) 



of a given age and metallicity (Worthey 1992 



1994). 



Using such models, Worthey showed that the age-metallicity 
degeneracy was actually worse than suspected: a factor of two 
uncertainty in the metallicity of a galaxy mimics a factor of 
three uncertainty in its age at fixed color or metal-line strength, 
the so-called "3/2 law". The law implies that such commonly 
used "age" indicators as colors and metal-line strengths are by 
themselves useless (although they still are widely used). At the 
same time, the Worthey models als o provided a quantitative too l 
to break the degeneracy (see also Worthey & Ottaviani 1997). 
A Balmer index plotted versus a metal line (or color) yields 
a two-dimensional theoretical grid; the equivalent single-burst 
age and metallicity for a population can be read off from its 
location in this grid. Tests of the method on composite stellar 
populations will be demonstrated in Trager et al. (1999; Pa- 
per II), where it is shown that these single-stellar-population 
(SSP) equivalent parameters correspond approximately to the 
luminosity-weighted vector addition of populations in the index 
diagrams. A galaxy's age determined from its integrated spec- 
trum is thus quite sensitive to recent star formation, and hence 
to the epoch and strength of its last major dissipative merger or 
accretion event. 

While Worthey models validated use of the Balmer lines, 
they also showed that extremely accurate Balmer data would be 
needed. To our knowledge, the line-strength data of Gonzalez 
(1993) are still the only published data on a diversified sample 
of local E galaxies that are adequate for this purpose. Apply- 
ing early Worthey models to his data, Gonzalez found that blue, 
weak-lined ellipticals in his sample tended to have young ages, 
while red, strong-lined ellipticals had older ages. In contrast, 
the metallicity spread was fairly small, le ss than a few tent hs of 
a dex. This result seemed to imply (G93, Faber et al. 1995) that 
age was the major cause of the well known color/line-strength 
relati on in the G93 sample, not metallicity as in the classic pic- 
ture ( iBaum 19591; |McClure & van den Bergh 1968|; S pinrad & 
Taylor 1971; [Faber 1972| , 1973). The large age spread in G93 
galaxies was later confirmed by Trager (1997) and by Tantalo, 
Chiosi & Bressan (1998a; TCB98) using later stellar population 
models. 

Excellent line strengths have also been measured for E and 
SO galaxies in the Fornax cluster by Kuntschner & Davies 
(1998) and Kuntschner (1998). Fornax turns out to be the re- 
verse of the G93 sample in showing a larger spread in metal- 
licity than age; the dense cluster environment of Fornax may 
be the key difference. A goal of the present series of papers 
is to explore the relative importance of age versus metallicity 



in driving the color and line strength relations of ellipticals in 
different environments (see Paper II). 

Although the data of G93 shed hope on solving the age vs. 
metallicity problem, they brought another simmering problem 
to the fore, namely, non-solar abundance ratios. Enhancement 
of Mg relative to Fe had been suggested by O'Connell (1976) 
and Peletier (1989) and shown to be widespread in the Lick/IDS 
ellipticals by Worthey, Faber & Gonzalez (1992). However, 
the high-quality data of G93 offered a great improvement for 
hard-to-measure weak Fe lines, and, using them, Trager (1997) 
showed that metallicities deduced from Mg were indeed consid- 
erably higher than those deduced from Fe. Other elements such 
as Na, C, N, a nd possibly O are also probably enhanced in gi- 
ant ellipticals (Worthey 1998 ). Because the Worthey models are 
not designed for non-solar abundance ratios, applying them to 
different metal line features in elliptical spectra gives inconsis- 
tent ages and, especially, abundances. The progress promised 
by the G93 data thus suddenly came to a full stop. 

The present paper addresses the problem of non-solar abun- 
dance ratios in a rough but hopefully satisfactory way. On the 
one hand, the general effects of non-solar ratios on evolutionary 
isochrones are no w beginning to be understood (Salaris , Chieffi 
& Straniero 19 93; [Weiss, Peletier & Matteucci 19951; Salaris 
& Weiss 1998; Bressan, priv. comm.; see Tantalo, Chiosi 
& Bressan 1998a). Second, the responses of nearly all the 
Lick/IDS indices to non-solar element ratios have been mod- 
eled by Tripicco and Bell (1995; TB95). The latter prove cru- 
cial, and it is really these responses that open the way forward. 
Using both inputs, reasonable corrections to the W94 indices 
for non-solar ratios can be estimated for the first time. The 
corrected models are used here to derive three SSP-equivalent 
population parameters for each galaxy — age, mean metallicity, 
and mean element "enhancement ratio." Future papers will use 
these parameters to study stellar populations as a function of 
galaxy type, determine correlations among age, metallicity, en- 
hancement, and other variables, and measure radial population 
gradients. 

Ot her groups ([Weiss, Peletier & Matteucci \995\, G reggio 
1997; [Trager 1997 ) have also attempted to interpret G93 data in 
terms of non-solar abundance ratios, but their approaches were 
more ad hoc. Inferred metallicities and enhancements both tend 
to be larger than what we find here. The most similar analysis 
so far is by Tantalo, Chiosi & Bressan (1998a; TCB98), build- 
ing o n previous work by that group (B ressan, Chiosi, & Tantalo 
1996). However, these authors use different response functions 
from ours (and in fact do not correct the Fe index at all for non- 
solar ratios). Their results consequently differ, and a section is 
devoted to comparing our results to their work (Sec. |6.1[ ). 

We note briefly that Balmer-line equivalent widths might be 
spuriously contaminated by light from blue horizontal-branch 
(BH B) stars or blue straggler stars (BSS) (e.g., Burstein et al. 
1984; [Lee 1994[ ; [Faber et al. 1995[; [Trager 1997| ). These possi- 
bilities are discussed in Section [5.1[ To anticipate the conclu- 
sions, we believe that current data do not support the existence 
of large numbers of BHB and BSS stars in giant elliptical galax- 
ies, and we thus conclude that the SSP-equivalent ages derived 
here for both young and old ellipticals must be substantially 
correct. Likewise, reduction of Balmer indices by emission fill- 
in, though present, cannot change the derived ages very much. 
Thus, despite efforts, we have been unable to find any expla- 
nation for the wide range of Balmer line strengths in the G93 
galaxies other than a wide range of SSP-equivalent ages. This 
is our principal conclusion. 
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The outline of this paper is as follows: SectionB presents 
absorption-line data for the G93 galaxies. Section ||presents a 
brief description of the Worthey (1994) models; their extension 
to non-solar abundance ratios using the results of TB95; the 
final choice of elements for inclusion in the enhanced element 
group; the method for determining the stellar population param- 
eters from the models; and the final population parameters for 
the G93 sample. Section || briefly presents the parameters for 
the G93, both central and global, and their distributions. Sec- 
tion H discusses the assumptions, in particular the use of H/3 
as an age indicator, and examines all known uncertainties in 
the age, metallicity, and abundance-ratio scales and zeropoints. 
Section ^presents evidence from other absorption-line strength 
studies for the presence of intermediate-age stellar populations 
in elliptical galaxies; it also compares in detail our results to 
those of TCB98. Two appendices discuss the effect of chang- 
ing isochrones in the models and the effect of using different 
prescriptions for emission and velocity dispersion corrections 
to H/3. 

2. DATA 
2. 1 . The galaxy sample 

The G93 galaxy sample was not selected according to quan- 
titative criteria but was rather chosen with the aim of covering 
relatively uniformly the full range of color, line strength, and 
velocity dispersion shown by local elliptical galaxies. As such, 
it contains more dim, blue, weak-lined, low-dispersion galax- 
ies than would be found in a magnitude-limited sample. In that 
sense the G93 sample may more closely resemble a volume- 
limited sample, but this has not been established quantitatively. 

The original sample in G93 consisted of 41 galaxies, of 
which 40 are included here. NGC 4278 has been discarded 
because of its strong emission. Table |l] presents morphologies, 
positions, and heliocentric redshifts. All galaxies are classi- 
fied as elliptical (o r compact elliptical) in the RC3 (de V au- 
couleurs et al. 1991) , the RSA (jSandage & Tammann 1987| ), or 
the Carnegie Adas ( |Sandage & Bedke 1994| ) except for NGC 
507 and NGC 6703, both classified as SAO in the RC3 but not 
cataloged in the RSA or the Carnegie Atlas. NGC 224 (the 
bulge of M 31) is also included. 

The environmental distribution of the G93 sample bears com- 
ment. Group assignments and approximate group richnesses 
may be found for nearly all galaxies in Faber et al. (1989). Most 
of the G93 galaxies are in poor groups, a few are quite isolated 
(there are no other galaxies in the RC3 within 1 Mpc projected 
distance and ±2000 km s"' of NGC 6702, for example [Col- 
bert, Mulchaey & Zabludoff, in prep.]), and six are members 
of the Virgo cluster Only one is in a rich cluster (NGC 547, 
in Abell 194). We therefore refer to the galaxies in this sample 
as local "field" ellipticals, given the low-density environments 
of most of them. Environmental effects are discussed in more 
detail in Paper II. 

2.2. G93 indices: calibrations and corrections 

The Lick/IDS indices were introduced by Burstein et al. 
(1984) to measure prominent absorption features in the spectra 
of old stellar populations in the 4100-6300 A region. A large 
and homog eneous database of s t ellar and galaxy sp ectra was 
assembled (Worthey et al. 1994; Frager et al. 1998, hereafter 
TWFBG98) with the Image Dissector Sc anner at Lick Obser- 
vatory (IDS; Robinson & Wampler 1972 ). A description of the 
Lick/IDS system and its application to stellar and galaxy spec- 



tra is given in those papers. 

Gonzalez (1993) measured Lick/IDS indices with a differ- 
ent spectrograph setup, at higher dispersion, and over a re- 
stricted spectral range (4700-5500 A). The four best indices 
in his wavelength interval are H/3, Mgb, Fe5270, and Fe5335, 
which we use in this paper The bandpasses of these four in- 
dices are given in Table g, and the precise index definitions are 
given in G93, Worthey et al. (1994), and TWFBG98. 

We use a combined "iron" index, (Fe), in this work, which 
has smaller errors than either Fe index separately and is defined 
as follows: 

, , Fe5270 + Fe5335 

(Fe) = . (1) 

It has the conve nient property of being sensitive primarily to 
[Fe/H] (see Sec. 3.1.2 ). Although Mg2 has also become a stan- 
dard "metallicity" indicator for the integrated spectra of galax- 
ies, we do not use it to determine stellar population parameters. 
G93 was unable to transform his observations of this broad in- 
dex (or of Mgi) accurately onto the Lick/IDS system due to 
chromatic focus variations in his spectrograph, coupled with 
the steep light gradient in the central regions of most ellipticals 
(Fisher et al. 1995 avoided Mg2 for the same reason). We pre- 
fer to use the narrower index Mgb, which is not affected by this 
problem. 

2.2.1. Velocity-dispersion corrections 

The observed spectrum of a galaxy is a convolution of the in- 
tegrated spectrum of its stellar population with the line-of-sight 
velocity distribution function of its stars. Indices measured 
for broad-line galaxies are therefore too weak compared to un- 
broadened standard stars. TWFBG98 statistically corrected the 
Lick/IDS indices for this effect in the following way: individ- 
ual stellar spectra of a variety of spectral types (plus M 32) 
were convolved with Gaussian broadening functions of increas- 
ing widths and their indices were remeasured. A smooth mul- 
tiplicative correction as a function of velocity dispersion was 
determined separately for each index and applied to the galaxy 
data. 

G93 used a more sophisticated technique, taking advantage 
of the higher resolution and signal-to-noise of his data. His 
stellar library was used to synthesize a summed stellar template 
representing a best fit to the the spectrum of each galaxy. In- 
dices were measured from the unbroadened template and again 
from the broadened template, generating a velocity dispersion 
correction for each galaxy that was tuned to its spectral type. 
For Mgb, Fe5270, and Fe5335, the mean multiplicative correc- 
tions of G93 are very similar to those of TWFBG98 (compare 
his Figure 4.1 with Figure 3 of TWFBG98). However, for H/3, 
the correction of G93 is flat or even negative, whereas the cor- 
rection of TWFBG98 is always positive and reaches the value 
1.07 at cr = 300 km s"'. Use of the TWFBG98 correction in- 
creases H/3 over G93 and leads to slightly younger ages. In 
what follows, we use the G93 correction to remain consistent 
with his published data but explore the effects of the TWFBG98 
correction in Appendix ^. The data marginally appear to favor 
TWFBG98, but the differences are not large. 

2.2.2. Emission corrections 

G93 noted that [O III] AA4959,5007 are clearly detectable 
in about half of the nuclei in his sample and that most of 
these galaxies also have detectable H/3 emission (see his Fig- 
ure 4.10). For galaxies in his sample with strong emis- 
sion, H/3 is fairly tightly correlated with [O III] such that 
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Table 1 

The Gonzalez (1993) Sample' 





Type 


Type 






cz 


Other 


Name 


(RC3) 


(CA) 


0^2000 


(^2000 


(km s-i) 


names 


NGC 221 


cE2 


cE2 


00:42:41.9 


+40:51:52 


-204 ± 7 


M32 


NGC 224 


SA(s)b 


Sb 


00:42:44.2 


+41:16:08 


-300 ± 7 


M31 


NGC 315 


E+: 




00:57:48.9 


+30:21:09 


4942 ± 6 




NGC 507 


SA(r)0 




01:23:39.8 


+33:15:23 


4908 ± 1 1 




NGC 547 


El 




01:26:00.7 


-01:20:44 


5468 ± 6 




NGC 584 


E4 


S0i(3,5) 


01:31:20.7 


-06:52:06 


1866 ± 6 




NGC 636 


E3 


El 


01:39:06.5 


-07:30:46 


1860 ± 6 




NGC 720 


E5 


E5 


01:53:00.4 


-13:44:18 


1741 ±11 




NGC 821 


E6? 


E6 


02:08:21.0 


+10:59:44 


1730 ± 7 




NGC 1453 


E2-3 


EO 


03:46:27.2 


-03:58:09 


3886 ± 6 




NGC 1600 


E3 


E4 


04:31:39.9 


-05:05:10 


4688 ± 8 




NGC 1700 


E4 


E3 


04:56:56.3 


-04:51:52 


3895 ± 7 




NGC 2300 


SAO 


E3 


07:32:22.0 


+85:42:27 


1938 ± 7 




NGC 2778 


E 




09:12:24.4 


+35:01:38 


2060 ± 7 




NGC 3377 


E5-6 


E6 


10:47:41.6 


+13:59:00 


724 ± 7 




NGC 3379 


El 


EO 


10:47:49.5 


+12:34:57 


945 ± 7 


M105 


NGC 3608 


E2 


El 


11:16:58.7 


+18:08:57 


1222 ± 7 




NGC 3818 


E5 


E5 


11:41:57.5 


-06:09:21 


1708 ±10 




NGC 4261 


E2-3 


E3 


12:19:23.2 


+05:49:31 


2238 ± 7 




NGC 4374 


El 


El 


12:25:03.7 


+12:53:14 


1060 ± 6 


M84 


NGC 4472 


E2 


El/SOi(l) 


12:29:46.5 


+07:59:48 


980 ±10 


M49 


NGC 4478 


E2 


E2 


12:30:17.4 


+12:19:44 


1365 ± 7 




NGC 4489 


E 




12:30:52.2 


+16:45:31 


970 ±10 




NGC 4552 


E 


SOi(O) 


12:35:39.9 


+12:33:25 


364 ± 7 


M89 


NGC 4649 


E2 


S0i(2) 


12:43:39.7 


+11:33:09 


1117± 6 


M60 


NGC 4697 


E6 


E6 


12:48:35.8 


-05:48:00 


1307 ±10 




NGC 5638 


El 


El 


14:29:40.4 


+03:14:04 


1649 ± 6 




NGC 5812 


EO 


EO 


15:00:57.0 


-07:27:19 


1929 ± 7 




NGC 5813 


El-2 


El 


15:01:11.2 


+01:42:08 


1954 ± 7 




NGC 5831 


E3 


El 


15:04:07.2 


+01:13:15 


1655 ± 5 




NGC 5846 


EO-1 


E4/S0i(4) 


15:06:29.3 


+01:36:21 


1714± 5 




NGC 6127 


E 




16:19:11.9 


+57:59:03 


4700 ±10 


NGC 6125/6128 


NGC 6702 


E: 




18:46:57.6 


+45:42:20 


4728 ± 5 




NGC 6703 


SAO- 




18:47:18.9 


+45:33:02 


2403 ± 7 




NGC 7052 


E 




21:18:32.9 


+26:26:48 


4672 ± 8 


B2 21 16+26 


NGC 7454 


E4 




23:01:06.6 


+16:23:24 


2051 ± 7 




NGC 7562 


E2-3 




23:15:57.4 


+06:41:16 


3608 ± 5 




NGC 7619 


E 


E3 


23:20:14.7 


+08:12:23 


3762 ± 5 




NGC 7626 


E pec: 


El 


23:20:42.4 


+08:13:02 


3405 ± 4 




NGC 7785 


E5-6 


S0i(5)/E5 


23:55:19.1 


+05:54:53 


3808 ± 5 





^NGC 4278 has been removed because of its strong emission lines. 



Note. — Col. (1): Galaxy name. Col. (2): Morphological type from RC3. Col. (3): Morpho- 
logical type from Carnegie Atlas (Sandage & Bedke 1994) or RSA (Sandage & Tammann 1987). 
Cols. (4)-(5): Position (J2000.0) from NED. Col. (6): Heliocentric radial velocity from Gonzalez 
(1993). Col. (7): Other galaxy names. 
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Table 2 

Lick/ids Indices Used in This Study 



./■ 


Name 


Index Bandpass 


Pseudocontinua 


Units 


Measures" 


09 


H/3 


4847.875^876.625 


4827.875-4847.875 


A 


H/3,(Mg) 








4876.625-4891.625 






13 


Mgb 


5160.125-5192.625 


5142.625-5161.375 


A 


Mg,(C),(Cr),(Fe) 








5191.375-5206.375 






14 


Fe5270 


5245.650-5285.650 


5233.150-5248.150 


A 


Fe,C,(Mg) 








5285.650-5318.150 






15 


Fe5335 


5312.125-5352.125 


5304.625-5315.875 


A 


Fe,(C),(Mg),Cr 



5353.375-5363.375 



"Dominant species; species in parentheses control index in a negative sense (index weak- 
ens as abundance grows). See Tripicco & Bell (1995) and Worthey (1998). 



EW(H/3em)/EW([0 III]A5007) - 0.7. A statistical correction 
of 



AH/3 = 0.7 EW([0 III]A5007) 



(2) 



was therefore added to H/? to correct for this residual emission. 

We have examined the accuracy of this correction by study- 
ing H/3/[0 III] among the G93 galaxies, supplemented by addi- 
tional early-type galaxies from the emission-line catalog of Ho, 
Filipenko & Sargent (1997). The sample was restricted to in- 
clude only normal, non-AGN Hubble types E through SO-, and 
to well measured objects with EW(Ha) > 1 .0 A. For 27 galax- 
ies meeting these criteria, H/3/[0 III] varies from 0.33 to 1.25, 
with a median value of 0.60. This suggests that a better correc- 
tion coefficient in Equation ^ might be 0.6 rather than 0.7, and 
thus that the average galaxy in G93 is slightly overcorrected. 
For a median [O III] strength through the G93 re/8 aperture of 
0.17 A, the error would be about 0.02 A, or 3% in age. This 
systematic error for a typical galaxy is negligible compared to 
other sources of error in the ages (see Table Random errors 
due to scatter in the ratio are about three times larger but are 
still small. 

Carrasco et al. (1996) report no correlation between H/3 and 
[O III] emission in their sample of early-type galaxies, but give 
no data. Their claim is explored in Appendix B, which repeats 
our calculations but with no H/3 correction. The ages of a few 
strong- [O III] galaxies are increased, as expected, but the broad 
conclusions of this work are unaffected. 

No correction for [N II] emission has been made to Mgb, al- 
though this has been suggested as a sometimes significant con- 
tributor to this index (by increa sing the flux in the red sideband; 
Goudfrooij & Emsellem 1996] ). Only NGC 315 and NGC 1453 
would be affected (see G93). 

Table |3^ presents final corrected index strengths, velocity 
dispersion corrections, and emission corrections for measure- 
ments through a central r(./8 aperture; Table |3^ presents simi- 
lar data for a global r^/l aperture. All values are taken directly 
from G93. The aperture index strengths are weighted averages 
of the major and minor axis profile data, computed so as to 
mimic what would be observed through the indicated circular 
aperture (see G93 for details). 

-There is a systematic color offset in the Kurucz (1992) models when compared 
red by 0.06 mag in (B-V) (but not in other colors; W94). All model (B-V) colors 



3. SSP-EQUIVALENT STELLAR POPULATION PARAMETERS 

3.1. Method 

3.1.1. Solar-abundance ratio models of Worthey (1994) 

SSP-equivalent population parameters have been derived 
by matching observed line strengths of Mgb, (Fe), and H/3 
to updated single-b urst stellar population (SSP) models of 
W94 (available at ^ttp://astro.sau.edu/^worthey/ ; "Padova" 
isochrones by Bertelli et al. (1994) are explored in Ap- 
pendix The models of W94 depend on two adjustable 
parameters — metallicity and single-burst age — and one fixed 
parameter, the initial mass function exponent (IMF), here cho- 
sen to have the Salpeter value. For reasons stated below, we be- 
lieve that the basic models of W94 have essentially solar abun- 
dance ratios; we will presently adjust these models to allow for 
non-solar ratios and, in the process, derive a third adjustable 
parameter, the non-solar enhancement ratio, [E/Fe]. The W94 
models are reviewed briefly here, and the reader is referred to 
Worthey (1994) for more details. 

The models incorporate three ingredients: stellar evolu- 
tionary isochrones, a stellar SED library, and absorption-Une 
strengths. From the bottom of the main sequence to the base 
of the red-giant branch (RGB), th e models use the is ochrones 
of VandenBerg a nd collaborators ([VandenBerg 1985]; Vanden- 
Berg & Bell 1985; [VandenBerg & Laskarides 1987p rhese are 
mated to red giant branches from th e Revised Yale Isochrones 
( Green, Demarque & King 1987 ) by shifting the latter in 
AlogL and A log to match at the base of the RGB. Extrapo- 
lations are made to cover a wide range of (Z, Y, age) assuming 
that Zq = 0.0169 and Y = 0.228 + 2.7Z. 

The SED library was constructed using the model atmo- 
spheres and SEDs of Kurucz (1992) for stars hotter than 3750 
K, and model SEDs of Bessel et al. (1989, 1991) and observed 
SEDs from Gunn & Stryker (1983) for cooler M giants.^ 

Polynomial fitting functions from Worthey et al. (1994) 
for the Lick/IDS indices are used as the basis of the model 
absorption-line strengths. Metal-rich stars in the Lick/IDS li- 
brary are a random sample of metal-rich stars in the solar neigh- 
borhood; since evidence suggests that such stars have essen- 

with the empirical colors of Johnson (1966), the Kurucz (1992) models being too 
in this series are corrected for this offset. 
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Table 3 a 

Fully Corrected Index Strengths in the Central r<,/8 Aperture' 



Name 


LO IIIJ 


(7 


Hp 


(T 


JVlg o 


(7 


rCDZ 


(7 




(7 


(Fe) 


(7 


NGC 221 


0.11 


0.05 


2.31 


0.05 


2.96 


0.03 


2.88 


0.04 


2.61 


0.04 


2.75 


0.03 


NGC 224 


0.11 


0.07 


1.67 


0.07 


4.85 


0.05 


2.88 


0.04 


2.61 


0.04 


3.09 


0.04 


NGC 315 


0.45 


0.06 


1.74 


0.06 


4.84 


0.05 


2.92 


0.06 


2.85 


0.07 


2.88 


0.05 


NGC 507 


0.05 


0.09 


1.73 


0.09 


4.52 


0.11 


2.95 


0.12 


2.60 


0.15 


2.78 


0.09 


NGC 547 


0.39 


0.06 


1.58 


0.07 


5.02 


0.05 


2.97 


0.07 


2.66 


0.08 


2.81 


0.05 


NGC 584 


0.27 


0.05 


2.08 


0.05 


4.33 


0.04 


3.03 


0.04 


2.77 


0.04 


2.90 


0.03 


NGC 636 


0.07 


0.06 


1.89 


0.04 


4.20 


0.04 


3.19 


0.05 


2.87 


0.05 


3.03 


0.04 


NGC 720 


0.14 


0.09 


1.77 


0.12 


5.17 


0.11 


2.94 


0.12 


2.80 


0.14 


2.87 


0.09 


NGC 821 


0.03 


0.05 


1.66 


0.04 


4.53 


0.04 


3.08 


0.05 


2.81 


0.05 


2.94 


0.04 


NGC 1453 


0.89 


0.06 


1.60 


0.06 


4.95 


0.05 


2.96 


0.06 


2.99 


0.07 


2.98 


0.05 


NGC 1600 


0.15 


0.06 


1.55 


0.07 


5.13 


0.06 


3.01 


0.07 


3.10 


0.09 


3.05 


0.06 


NGC 1700 


0.22 


0.05 


2.11 


0.05 


4.15 


0.04 


3.17 


0.05 


2.83 


0.05 


3.00 


0.04 


NGC 2300 


0.08 


0.06 


1.68 


0.06 


4.98 


0.05 


3.04 


0.06 


2.89 


0.07 


2.97 


0.05 


NGC 2778 


0.69 


0.07 


1.77 


0.08 


4.70 


0.06 


3.01 


0.07 


2.69 


0.08 


2.85 


0.05 


NGC 3377 


0.39 


0.05 


2.09 


0.05 


3.99 


0.03 


2.77 


0.04 


2.44 


0.04 


2.61 


0.03 


NGC 3379 


0.18 


0.05 


1.62 


0.05 


4.78 


0.03 


2.98 


0.04 


2.73 


0.04 


2.85 


0.03 


NGC 3608 


0.18 


0.06 


1.69 


0.06 


4.61 


0.04 


3.13 


0.05 


2.75 


0.06 


2.94 


0.04 


NGC 3818 


0.33 


0.08 


1.71 


0.08 


4.88 


0.07 


3.09 


0.08 


2.85 


0.08 


2.97 


0.06 


NGC 4261 


0.21 


0.06 


1.34 


0.06 


5.11 


0.04 


3.14 


0.05 


2.88 


0.06 


3.01 


0.04 


NGC 4374 


0.37 


0.05 


1.51 


0.04 


4.78 


0.03 


2.94 


0.04 


2.69 


0.04 


2.82 


0.03 


NGC 4472 


0.01 


0.08 


1.62 


0.06 


4.85 


0.06 


2.97 


0.07 


2.84 


0.08 


2.90 


0.05 


NGC 4478 


0.06 


0.06 


1.84 


0.06 


4.33 


0.05 


3.03 


0.06 


2.84 


0.06 


2.93 


0.04 


NGC 4489 


0.11 


0.08 


2.39 


0.07 


3.21 


0.06 


2.89 


0.07 


2.44 


0.07 


2.67 


0.05 


NGC 4552 


0.25 


0.05 


1.47 


0.05 


5.15 


0.03 


3.02 


0.04 


2.95 


0.04 


2.98 


0.03 


NGC 4649 


0.09 


0.05 


1.40 


0.05 


5.33 


0.04 


3.01 


0.04 


3.01 


0.05 


3.01 


0.03 


NGC 4697 


0.10 


0.07 


1.75 


0.07 


4.08 


0.05 


2.97 


0.06 


2.57 


0.06 


2.77 


0.04 


NGC 5638 


0.00 


0.06 


1.65 


0.04 


4.64 


0.04 


3.02 


0.05 


2.66 


0.05 


2.84 


0.04 


NGC 5812 


0.05 


0.06 


1.70 


0.04 


4.81 


0.04 


3.09 


0.05 


3.02 


0.06 


3.06 


0.04 


NGC 5813 


0.27 


0.06 


1.42 


0.07 


4.65 


0.05 


2.83 


0.06 


2.52 


0.07 


2.67 


0.04 


NGC 5831 


0.18 


0.05 


2.00 


0.05 


4.38 


0.04 


3.17 


0.04 


2.92 


0.04 


3.05 


0.03 


NGC 5846 


0.39 


0.08 


1.45 


0.07 


4.93 


0.05 


2.95 


0.06 


2.77 


0.06 


2.86 


0.04 


NGC 6127 


0.04 


0.08 


1.50 


0.05 


4.96 


0.06 


2.90 


0.07 


2.79 


0.08 


2.85 


0.06 


NGC 6702 


0.40 


0.06 


2.46 


0.06 


3.80 


0.04 


3.02 


0.05 


2.97 


0.06 


2.99 


0.04 


NGC 6703 


0.36 


0.05 


1.88 


0.06 


4.30 


0.04 


3.06 


0.05 


2.79 


0.05 


2.92 


0.04 


NGC 7052 


0.43 


0.06 


1.48 


0.07 


5.02 


0.06 


2.89 


0.07 


2.78 


0.08 


2.83 


0.05 


NGC 7454 


0.11 


0.06 


2.15 


0.06 


3.27 


0.05 


2.68 


0.06 


2.27 


0.06 


2.47 


0.04 


NGC 7562 


0.09 


0.05 


1.69 


0.05 


4.54 


0.04 


3.08 


0.05 


2.65 


0.05 


2.87 


0.03 


NGC 7619 


-0.02 


0.05 


1.36 


0.04 


5.06 


0.04 


3.03 


0.05 


3.08 


0.06 


3.06 


0.04 


NGC 7626 


0.11 


0.05 


1.46 


0.05 


5.05 


0.04 


2.85 


0.05 


2.80 


0.05 


2.83 


0.03 


NGC 7785 


0.15 


0.06 


1.63 


0.06 


4.60 


0.04 


2.88 


0.05 


2.94 


0.06 


2.91 


0.04 



^From Gonzalez (1993), Table 4.7. 
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Table 3b 

Fully Corrected Index Strengths in the rgjl Aperture^ 



Name 


[0 111] 


a 


H/3 


(7 






Fe52 


(T 


Fe53 


a 


(Fe) 


a 


NGC 221 


0.14 


0.07 


2.15 


0.07 


2.96 


0.07 


2.79 


0.06 


2.47 


0.06 


2.63 


0.04 


NGC 315 


0.25 


0.07 


1.80 


0.08 


4.52 


0.09 


2.62 


0.09 


2.72 


0.11 


2.67 


0.07 


NCC 507 


0.27 


0.13 


2.06 


0.17 


4.69 


0.17 


2.72 


0.17 


2.18 


0.19 


2.45 


0.13 


NCC 547 


0.19 


0.08 


1.42 


0.10 


4.80 


0.09 


2.69 


0.09 


2.60 


0.11 


2.65 


0.07 


NCC 584 


0.19 


0.07 


2.06 


0.07 


4.13 


0.07 


2.77 


0.06 


2.55 


0.06 


2.66 


0.04 


NCC 636 


0.07 


0.07 


1.87 


0.06 


3.98 


0.08 


2.93 


0.07 


2.56 


0.07 


2.75 


0.05 


NCC 720 


0.32 


0.11 


2.28 


0.16 


4.98 


0.16 


2.78 


0.16 


2.92 


0.18 


2.85 


0.12 


NCC 821 


0.04 


0.07 


1.82 


0.06 


4.11 


0.08 


2.83 


0.08 


2.68 


0.08 


2.75 


0.06 


NCC 1453 


0.78 


0.07 


1.69 


0.09 


4.43 


0.08 


2.84 


0.08 


2.81 


0.09 


2.83 


0.06 


NCC 1600 


0.15 


0.08 


1.74 


0.09 


5.21 


0.10 


3.03 


0.10 


3.07 


0.12 


3.05 


0.08 


NCC 1700 


0.20 


0.07 


2.11 


0.07 


3.90 


0.08 


2.94 


0.07 


2.71 


0.07 


2.83 


0.05 


NCC 2300 


0.05 


0.07 


1.63 


0.06 


4.70 


0.09 


2.78 


0.08 


2.72 


0.09 


2.75 


0.06 


NCC 2778 


0.60 


0.08 


1.56 


0.10 


4.44 


0.09 


2.88 


0.09 


2.40 


0.10 


2.64 


0.07 


NCC 3377 


0.50 


0.07 


2.13 


0.07 


3.46 


0.07 


2.51 


0.06 


2.06 


0.06 


2.29 


0.04 


NCC 3379 


0.11 


0.07 


1.59 


0.06 


4.44 


0.07 


2.80 


0.06 


2.55 


0.06 


2.67 


0.04 


NCC 3608 


0.11 


0.07 


1.73 


0.09 


4.04 


0.08 


3.09 


0.07 


2.61 


0.08 


2.85 


0.05 


NCC 3818 


0.42 


0.10 


1.81 


0.11 


4.19 


0.12 


2.73 


0.11 


2.54 


0.11 


2.63 


0.08 


NCC 4261 


-0.03 


0.07 


1.30 


0.06 


4.75 


0.08 


3.05 


0.07 


2.53 


0.08 


2.79 


0.05 


NCC 4374 


0.22 


0.07 


1.56 


0.06 


4.50 


0.07 


2.78 


0.06 


2.57 


0.06 


2.67 


0.04 


NCC 4472 


0.03 


0.10 


1.67 


0.09 


4.60 


0.11 


2.83 


0.10 


2.80 


0.11 


2.81 


0.07 


NCC 4478 


0.04 


0.07 


1.73 


0.06 


4.17 


0.08 


2.74 


0.07 


2.52 


0.07 


2.63 


0.05 


NCC 4489 


0.03 


0.10 


2.27 


0.08 


2.83 


0.12 


2.77 


0.11 


2.14 


0.11 


2.46 


0.08 


NCC 4552 


0.15 


0.07 


1.52 


0.06 


4.80 


0.07 


2.83 


0.06 


2.70 


0.06 


2.77 


0.04 


NCC 4649 


0.15 


0.07 


1.38 


0.07 


5.13 


0.07 


2.62 


0.07 


2.75 


0.07 


2.69 


0.05 


NCC 4697 


0.01 


0.10 


1.66 


0.07 


3.62 


0.10 


2.55 


0.09 


2.28 


0.09 


2.42 


0.06 


NCC 5638 


-0.02 


0.07 


1.68 


0.06 


4.12 


0.08 


2.79 


0.07 


2.50 


0.07 


2.65 


0.05 


NCC 5812 


0.08 


0.07 


1.71 


0.05 


4.54 


0.08 


2.96 


0.07 


2.82 


0.07 


2.89 


0.05 


NCC 5813 


0.25 


0.08 


1.23 


0.10 


4.26 


0.09 


2.90 


0.09 


2.43 


0.10 


2.67 


0.07 


NCC 5831 


0.18 


0.07 


2.03 


0.07 


3.88 


0.07 


2.89 


0.06 


2.55 


0.06 


2.72 


0.04 


NCC 5846 


0.11 


0.10 


1.27 


0.10 


4.55 


0.11 


2.63 


0.10 


2.67 


0.10 


2.65 


0.07 


NCC 6127 


0.05 


0.10 


1.50 


0.07 


4.67 


0.11 


2.76 


0.10 


2.55 


0.11 


2.65 


0.07 


NCC 6702 


0.46 


0.07 


2.49 


0.09 


3.71 


0.08 


2.91 


0.08 


2.84 


0.09 


2.88 


0.06 


NCC 6703 


0.26 


0.07 


1.83 


0.07 


4.00 


0.08 


2.82 


0.07 


2.53 


0.07 


2.67 


0.05 


NCC 7052 


0.24 


0.08 


1.77 


0.09 


4.66 


0.09 


2.78 


0.09 


2.70 


0.10 


2.74 


0.07 


NCC 7454 


0.05 


0.07 


2.08 


0.06 


2.88 


0.09 


2.39 


0.08 


2.10 


0.08 


2.25 


0.06 


NCC 7562 


0.13 


0.07 


1.72 


0.07 


4.42 


0.08 


2.92 


0.07 


2.54 


0.07 


2.73 


0.05 


NCC 7619 


-0.03 


0.07 


1.47 


0.05 


4.70 


0.08 


2.78 


0.07 


2.84 


0.08 


2.81 


0.05 


NCC 7626 


-0.03 


0.07 


1.44 


0.06 


4.64 


0.08 


2.70 


0.07 


2.53 


0.08 


2.62 


0.05 


NCC 7785 


0.05 


0.07 


1.52 


0.06 


4.30 


0.08 


2.81 


0.07 


2.82 


0.08 


2.81 


0.05 



''From Conzalez (1993), Table 6.1. 
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tially solar ratios of O, Mg, Na, an d other key elements rela- 
tive to Fe ( Edvardsson et al. 1993), we assume that the line- 



Following previous practice, we adopt the convention that 
a certain group of elements is "enhanced" in elliptical galax- 
ies (more is said on this below). Precisely which elements 
are enhanced, and by how much, is poorly known. From an 
intercomparison of absorption-line strengths in the Lick/IDS 
galaxy sample (TWFBG98), Worthey (1998) suggested that 
Mg, Na, and N are enhanced in giant ellipticals b ut that Ca 
ti-acks Fe (cf. [O'Connell 197^ ; [Vazdekis et al. 1996| ). Compar- 
ing to additional galaxy data from TWFBG98 below, we sug- 
gest that C also belongs to the enhanced group. Unfortunately, 
the Lick/IDS system has no indices th at are capable o f directly 
probing oxygen in elliptical galaxies ( Worthey 1998| ). Oxygen 
is important because it dominates [Z/H] on account of its high 
masss fraction. 

Because O (and perhaps C) are uncertain, we have con- 
sidered four models for the enhancement pattern in ellipti- 
cal galaxies, as described in Table 0. In each model there 
are three groups of elements: enhanced, depressed, and fixed. 
The assignment of elements to the three groups is always the 
same except for C and O, whose assignments vary. Elements 
i n the fixed group have their sola r (photospheric) abundances 
(Grevesse, Noels & Sauval 1996), while elements in each of 



strengths produced by the metal-rich models of W94 reflect 
solar-abundance ratios. 

To construct a model of a given age and metallicity, the ap- 
propriate stellar isochrone is first selected. Each star on the 
isochrone is assigned an SED from the flux library and a set 
of absorption-line strengths from the Lick/IDS fitting func- 
tions. Final model outputs are the integrated fluxed SED (from 
which colors and magnitudes can be derived) and absorption- 
line strengths on the Lick/IDS system. 

The ability of the W94 models to break the age-metallicity 
degeneracy is illustrated in Figures |l] and^, which plot Mgb 
and (Fe) versus H/3 for the G93 galaxies. Model grids from 
W94 are overplotted. Both line strength pairs break the de- 
generacy, but (Fe)-H/3 does so more than Mgfe-H/3 because 
(Fe) is less temperature sensitive than Mgb. Metallicities in- 
ferred from Mgb are clearly higher than from (Fe), reflecting 
the probable element enhancement in [Mg/Fe]. Low-H/3 galax- 
ies tend to fall off the grid to high ages in the (Fe)-H/3 diagram, 
especially through the rg/2 aperture (Figure 2b). Most of this 
effect is removed when (Fe) is corrected for depressed [Fe/H] 
(see below), and any small remainder can be attributed to use 
of the G93 H/3 velocity corrections (instead of TWFBG98) and, 
in a few galaxies, to possible residual, uncorrected H/3 emission 
(see Appendix ^). 

3. 1 .2. Non-solar abundance ratio models 

Adjusting the W94 models for non-solar ratios involves two 
steps. First, one must compute new evolutionary tracks in a 
fully self-consistent manner using new interior opacities, re- 
action rates, and atmospheric boundary conditions that faith- 
fully reflect the altered compositions. Second, one must com- 
pute new absorption-line indices. Part one is less developed in 
the literature but also proves to be less important; we discuss 
it below. Part two, the indices, could in principle be handled 
by observing populations of stars with known non-solar ratios 
and deriving empirical fitting functions for them. For exam- 
ple, Borges et al. (1995) derived a fitting function for Mg2 ver- 
sus [Mg/Fe] using local dwarf and subgiant stars (this was the 
function adopted by TCB98 for their population models); and 
Weiss, Peletier & Mateucci (1995) attempted to correct Mg2 
and (Fe) using Galactic Bulge stars studied by Rich (1988). 

However, it is hard to identify groups of stars with exactly 
the same (known) enhancements, and it is even more diffi- 
cult to vary the pattern of element abundance enhancements 
in a controlled way using real stars. For these reasons, a 
theoretical approach is recommended, and we have chosen to 
utilize the computations of Tripicco & Bell (1995), who re- 
computed all of the Lick/IDS spectral indices from a grid of 
theoretical stellar SEDs and atmospheres with varying abun- 
dance ratios. For three sample locations on an old stellar 
isochrone, TB95 tabulate the response of each Lick/IDS index 
to separate enhancements [X/H] = +0.3 dex for the elements 
X = C,N,0,Mg,Fe,Ca,Na,Si,Cr and Ti. These response func- 
tions are the basis for our corrections to the indices for non- 
solar abundance ratios. Note that, because we use the response 
functions differentially, we are insensitive to any zeropoint un- 
certainties that the TB95 indices may have (which are in any 
case known to be small, as TB95 showed by comparing to real 
stars). 

^In retrospect it would have made more sense to group N with C since they are nucleosynthetically related (Woosley & Weaver 1995), but making this change 
would negligibly affect the conclusions. 



the enhanced and depressed groups are all scaled up or down 
by the same factor After the amount of enhancement is chosen 
and C and O are assigned to their proper groups, the depres- 
sion of the depressed elements is calculated so as to preserve 
constant [Z/H]. 

In the present work, we generally take the enhanced group 
to include the abundant elements that are nucleosynthetically 
related to Mg, several of which are actually seen to be over- 
abundant in giant elliptical galaxies (Worthey 1998). Elements 
placed in the enhanced group include N, Ne, Na, Mg, Si, and 

5 (plus sometimes C and/or O).^ The iron-peak elements Cr, 
Mn, Fe, Co, Ni, Cu and Zn constitute the depressed group. All 
other elements (including those heavier than Zn) are in the fixed 
group, with the exception of Ca (in the depressed group), and 
C and O (which vary). 

As noted, the four models differ in their treatment of C and 
O: model 1 has C fixed, O up; model 2 has C fixed, O fixed; 
model 3 has C down, O up; and model 4 has C up, O up. Be- 
cause O is produced in massive stars like Mg, it is probable that 
it, too, is routinely enhanced in giant ellipticals; hence model 2 
is unlikely on nucleosynthetic grounds. Model 3, with C down 
and O up, is similar to the models of Weiss, Peletier & Mateucci 
(1995), TCB98, and Salaris & Weiss (1998). We show below 
that depression of C does not match the Lick/IDS indices and 
that this model is also therefore unlikely. Model 4, with C and 
O both enhanced, is our preferred model based on Mc William 

6 Rich (1990) and Rich & McWilHam (priv. comm.), who find 
that O and C are enhanced in lockstep with Mg in stars in the 
Galactic Bulge. However model 1 (with C fixed) is very hard 
to distinguish observationally from model 4 (see below). 

Because the enhanced elements are not exactly the same as 
the a-elements (e.g., Ca is nominally an a-element but appar- 
ently tracks Fe in elliptical spectra; Worthey 1998; TWFBG98), 
we use the notation [E/Fe], where "E" refers to the mass frac- 
tion of elements that are specifically enhanced in each model, 
in preference to the more common notation [a/Fe] used by pre- 
vious authors. Following TCB98, we write 



[Fe/H] = [Z/H] -A[E/Fe], 



(3) 
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0.5 0.6 
log Mgb 



0.5 
log (Fe) 



Fig. 1. — Line strengths of early-type galaxies in the Gonzalez (1993) sample in the central r(,/8 aperture. Model grids from Worthey (1994) are superimposed: 
solid lines are contours of constant age (from top, 1, 1.5, 2, 5, 8, 10, 15, 18 Gyr), and dotted lines are contours of constant [Z/H] (from left, [Z/H] = —0.5, —0.25, 0., 
0.25, 0.5, except at ages younger than 8 Gyr, where from left [Z/H] = -0.225, 0., 0.25, 0.5). (a) Mg^; versus H/3. (b) (Fe) versus H/3. Differences in the ages and 
metallicities inferred from these two diagrams result from the non-solar abundance ratios of giant elliptical galaxies. 




0.5 0.0 
log Mgb 



0.5 

log <Fe> 



Fig. 2. — Line strengths of early-type galaxies in the Gonzalez (1993) samp 
global loci here are shifted horizontally to the left of the central loci in Figure 



B through the global rejl aperture. The remainder of the figure is as in Figure |l|. The 
reflecting the fact that the outer parts are both older and more metal-poor. 
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Table 4 

Element Mass Fractions in Non-solar Abundance Ratio Models" 



Model 




A 


A/(l-A) 


A[E/Fe] 


A[E/H] 


/(C) 


fiO) 


/(Fe) 


/(E) 


[^hpe/^lpe] 


1 


CoO+ 


0.914 


10.6 


0.323 


0.023 


0.173 


0.509 


0.041 


0.267 


0.061 


2 


CoOo 


0.774 


3.42 


0.365 


0.065 


0.173 


0.482 


0.041 


0.294 


0.020 


3 


C-0+ 


0.759 


3.15 


0.369 


0.069 


0.087 


0.565 


0.041 


0.297 


0.001 


4 


C+0+ 


0.929 


11.7 


0.319 


0.019 


0.181 


0.504 


0.041 


0.265 


0.067 


Solar 












0.172 


0.482 


0.082 


0.253 


0.000 



"At fixed A[Fe/H] = -0.3 dex and [Z/H] = 0. 



Note. — Cols. (l)-(2): Enhancement model. All models have enhanced N, Ne, Na, Mg, Si, and S. Ca and 
the Fe-peak elements (Cr, Mn, Fe, Co, Ni, Cu, Zn) are depressed. All other elements are fixed at their so- 
lar abundances except that model 1 has C fixed at its solar abundance and O enhanced (CoO+); model 2 has 
both C and O fixed at their solar abundances (CoOo); model 3 has C depressed like the Fe group and O en- 
hanced (C-0+; this is the enhancement pattern favored by TCB98 and by Salaris & Weiss 1998); and model 4 
has both C and O enhanced (C+0+). "Solar" is taken from the solar photospheric abundances of Grevesse et 
al. (1996). Cols. (3)-(4): A is the response of the "enhanced" elements to changes in [Fe/H] at fixed [Z/H]: 
A[Fe/H] =-AA[E/Fe] =-A/(l-A)A[E/H]. Col. (5): A[E/Fe] at A[Fe/H] =-0.3 dex. Col. (6): A[E/H] at 
A[Fe/H] = -0.3 dex. Col. (7): Mass fraction of metals in C. Col. (8): Mass fraction of metals in O. Col. (9): 
Mass fraction of metals in Fe peak elements. Col. (10): Mass fraction of metals in enhanced elements (except for 
C and O). Note that Cols. (7)-(10) do not add up precisely to 1, as elements fixed at their solar abundances (other 
thanC and O) are not included. Col. (11): The ratio [Xhpe/^lpe] = [-'''c+-'''N+-^o+-^Ne/-^Mg+-'''si+-'''s+-'''ca+-'''Fe] 
(Salaris, Chieffi & Straniero 1993); see text. 



or 



A[Fe/H] = -AA[E/Fe]: 



1-A 



A[E/H] 



(4) 



at constant [Z/H], where their very small second-order term in 
[E/Fe] has been ignored. Table ^ gives values of A and illustra- 
tive heavy-element fractions (C, O, E-group, Fe-peak) for the 
four models, all at A[Fe/H] =-0.3 dex and solar [Z/H]; values 
of [E/Fe] and [E/H] for other values of A[Fe/H] can be cal- 
culated using Eqs. |] andQ For reference, TCB98's model has 

^TCB98 = 0.8. 

Table ^reveals an important fact — ^because the Fe-peak con- 
tribution to [Z/H] is so small (only 8% for solar abundance), 
reducing it by even 0.3 dex frees up only minimal room for the 
so-called "enhanced" elements. Hence, what really happens in 
enhanced models is that the enhanced (and fixed) elements re- 
main nearly at their solar values, whereas Fe ( and related ele- 
ments) are depressed. In short, we should think of giant ellipti- 
cals as failing to make Fe-peak elements rather than making too 
much of certain other elements. Likewise, the quantity [E/Fe] 
is not really an enhancement of the E-elements but rather a de- 
pression ofFe. 

Other autho rs, including ourselves (e.g., Worthey , Faber & 
Gon zalez 1992; [Weiss, Peletier & Matteucci 19951; G reggio 
1 997; |Vazdekis et al. 1997| ) have said this, but the contradic- 
tory notion nevertheless persists that strong Mg indices are due 
to an "overabundance" of Mg — this is not mathematically pos- 
sible if Mg, O, and the a-elements track one another closely, 
as these elements together dominate [Z/H] by mass. We show 
below that the TB95 response functions provide an alternative 
means of strengthening Mgb and Mg2, namely, via weak Fe- 
peak elements (see below). This unanticipated a«f /-correlation 



between Mgb and the Fe-peak elements is one of the major new 
features of our treatment and the cause of our relatively small 
derived values of [E/Fe] (compared with previous authors; see 
Sec. |). 

We return next to the problem of the stellar evolutionary 
isochrones. Since a full library of isochrones is not available 
for all abundance ratios, we follow the lead of TCB98, who 
suggest from examining their unpublished isochrones that mod- 
els with varying [E/Fe] are "virtually indistinguishable in the 
CMD" from models at the same [Z/H] with [E/Fe] = 0. Earlier, 
Salaris et al. (1993) had shown (at sub-solar metallicities) that 
a-enhanced isochrones are identical to scaled-solar abundance 
isochrones at the same Z provided that the quantity 



(5) 



^HPE 




+ '^N + ^0 + ^Ne 


. '^LPE . 




XMg + Xsi+Xs+ Xca + Xpe 



remains constant at the solar value (= 0). Here X, is the mass 
fraction in element /, and brackets indicate the usual logarithm 
relative to solar. The elements in Xhpe have high ionization 
potentials and their opacity governs the mean turnoff temper- 
ature; the elements in Xlpe have low ionization potentials and 
their opacity governs the temperature of the giant branch. Pre- 
serving the ratio [^hpe/^lpe] thus preserves the shape of the 
track, they say, and the new track is found to fit neatly into the 
old sequence at the same value of Z. Values of [Xhpe/^lpe] 
are given in Table || for our four models. Models 2 and 3 are 
nearly solar, while models 1 and 4 are about 15% overabundant 
in HPE elements. These small deviations prove to be relatively 
unimportant, as shown in Section 5.4. 

More recently, Salaris & Weiss (1998) have suggested that, at 
higher metallicities near solar, track constancy may break down 
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Table 5 

Fractional index responses to non-solar abundance ratios^ 



Model 


Star 


%F 


5C2 


5HI3 


<5Mgi 


(5Mg2 


SMgb 


(5Fe52 


(5Fe53 


<5(Fe) 


1 CoO+ 


eg 


53 


-0.037 


0.000 


0.128 


0.088 


0.295 


-0.113 


-0.191 


-0.147 




to 


44 


-0.130 


0.026 


0.617 


0.025 


0.124 


-0.226 


-0.260 


-0.240 




cd 


3 


-0.193 


-0.361 


0.183 


0.086 


0.153 


-0.177 


-0.210 


-0.192 




total 




-0.046 


0.027 


0.146 


0.080 


0.249 


-0.136 


-0.203 


-0.165 


2 CoOo 


eg 


53 


-0.043 


0.000 


0.179 


0.129 


0.366 


-0.124 


-0.205 


-0.159 




to 


44 


-0.176 


0.031 


0.541 


0.049 


0.184 


-0.226 


-0.261 


-0.241 




cd 


3 


-0.242 


-0.317 


0.238 


0.118 


0.192 


-0.191 


-0.228 


-0.208 




total 




-0.056 


0.031 


0.193 


0.118 


0.315 


-0.145 


-0.215 


-0.176 


3 C-0+ 


eg 


53 


-0.703 


0.000 


-0.342 


0.010 


0.767 


-0.198 


-0.138 


-0.172 




to 


44 


-0.727 


0.027 


-0.558 


-0.062 


0.296 


-0.232 


-0.281 


-0.252 




cd 


3 


-0.864 


0.314 


-0.021 


0.023 


0.168 


-0.199 


-0.208 


-0.203 




total 




-0.708 


0.026 


-0.330 


0.001 


0.625 


-0.204 


-0.165 


-0.187 


4 C+0+ 


eg 


53 


0.055 


0.000 


0.175 


0.095 


0.265 


-0.107 


-0.195 


-0.145 




to 


44 


-0.049 


0.026 


0.786 


0.032 


0.111 


-0.226 


-0.259 


-0.240 




cd 


3 


-0.075 


-0.395 


0.200 


0.091 


0.152 


-0.175 


-0.210 


-0.191 




total 




0.045 


0.027 


0.195 


0.086 


0.225 


-0.131 


-0.206 


-0.164 


TCB98 


total'' 






-0.050 




0.172 








-0.125 



"For a W94 model with age 12 Gyr, fixed A[Fe/H] = -0.3 dex, and [Z/H] = 0. 
''For the TCB98 model with an age of 10 Gyr and [Z/H] = 0. 

Note.— Cols. (l)-(2): Enhancement model from Table 6; "TCB98" is close to the model of Tantalo et al. (1998a; 
TCB98); model 4 is our prefeired model. Col. (3): Stellar component on isochi'one, following TB95. "eg" refers to cool 
giants (Te = 4255 K, logg = 1 .9); "to" refers to turn-off stars (Te = 6200 K, logg = 4.1); "cd" refers to cool dwarfs (Te = 4575 
K, logg = 4.6); and "total" refers to the composite spectrum. Col. (4): Typical fractional light contiibution from each stel- 
lar component to total flux, in percent (C2 has a slightly higher contribution from turn-off stars than the other indices). 
Col. (5)-(12): Fractional index responses for A[Fe/H] = —0.3 dex, in the sense SI = A///, where A/ is the index change 
and / is the original valuepf the index. These are computed by inserting the response functions of Tripicco & Bell (1995) 
for each element into Eq. p. (The thi'ee stellar entries do not sum to the total value because they are weighted by the unequal 
stellar indices of each type.) Fractional responses for the TCB98 model are based on element abundances computed from 
their Eqs. (3)-(5) and their Table 1. 
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and that increasing [a/Fe] both shifts the track to the blue and 
changes its shape. It is not clear whether these effects are due 
to high [E/Fe], to [Xhpe/^lpe] 7^0, or both . However, the mo- 
tions are small, and we show in Section 5.4 that their impact on 
the indices is probably slight. 

If isochrones do not shift (at fixed metallicity), we can as- 
sume that logg, logTf,, logL, and the SED of each star on the 
track are also constant. Hence, it is necessary only to calculate 
the changes in each spectral feature using the index response 
functions of TB95, by perturbing each element up or down ac- 
cording to the model. TB95 tabulate fractional index changes 
for three typical stars, one on the lower main sequence, one at 
the turnoff, and one on the RGB, at solar metallicity. We as- 
sume the same fractional changes at all metallicities and com- 
bine these responses by weighting by the fractional light con- 
tributions of each type of star at each index.** Details are given 
in the notes to Table^ 

Note that the TB95 response functions are for enhancement 
values corresponding to [X/H] = +0.3 dex. Response functions 
for arbitrary values of [E/Fe] are calculated via Eq. ^ to get 
A[Fe/H] and A[E/Fe] and then by exponentially scaling the 
response functions in Tables 4-6 of TB95 by the appropriate 
element abundance. The fractional response of index / is there- 
fore 



free parameters requires three indices, for which we use H/3, 
(Fe), and Mgb. First, an expanded model grid of line strengths 
as a function of t, [Z/H], and (now) [E/Fe] is generated by ap- 
plying the TB95 response functions to the base W94 models at 
each (f, [Z/H]). These new grids (one for each model 1-4) are 
created by interpolating the W94 models at intervals of Af = 0. 1 
Gyr and A[Z/H] = 0.01 and then interpolating the TB95 results 
at intervals of A[E/Fe] = 0.01 at each (f , [Z/H]). The process 
is then inverted to derive (f , [Z/H], [E/Fe]) for each galaxy by 
searching in the grid to find that point with minimum distance 
from the observed parameters {Hf3,Mgb, (Fe)). It was neces- 
sary to linearly extrapolate the W94 models to slightly higher 
ages and to both lower and higher [Z/H] values to cover the full 
range of (HP, Mgb, (Fe})-space populated by the observations. 
The range of (f , [Z/H], [E/Fe]) space covered by the final grids 
is 



1 < t (Gyr) < 22, 



22 < t (Gyr) < 30, 



-0.5 < [Z/H] < 1.25, 

-0.3 < [E/Fe] < 0.75 

-0.5 < [Z/H] < 0.5, 

-0.3 < [E/Fe] < 0.75. 



^=|n[l+^o.3(X,)] 



([Xi/H]/0.3) 



1, 



(6) 



where Rq 3(X,) is the TB95 response function for element / at 
[X//H]=+0.3dex.5 

Table || shows changes in the indices corresponding to the 
four models in Table all of which have A [E/Fe] ^0.3. H/3 is 
virtually unaffected by non-solar abundance ratios, even at sub- 
stantial [E/Fe]; all changes are less than 3%, which translates 
to < 8% in age. Changes in (Fe) are roughly the same in all 
models and amount to a decrease of about 20%, driven mostly 
by the decrease in [Fe/H] (of 0.3 dex). However, C24668, Mgi, 
Mg2, and Mg b are all different, owing to the presence (or not) 
of C2 bands in the passband or sidebands of these indices; C2 
and Mgi increase greatly with increasing C, Mgb declines with 
increasing C, while Mg2 stays about the same independent of 
C. These changes all reflect the different abundance of C in the 
models since the abundance of Mg (and other elements in the E 
group) is always about constant (cf. Table 

Finally, we note that Mg/? increases in all models, in apparent 
contradiction to the near constancy of [E/Fe]. This increase is 
due mostly to the decrease in Fe and Cr (see TB95), which has 
the effect of increasing Mgb. In fact, changes in all the Mg in- 
dices are driven more by the Fe-peak deficit than by any actual 
increase in Mg, proving once again that a more correct way of 
looking at elliptical galaxies is to regard them as Fe-poor rather 
than a-enhanced. 

3.2. Ages, metallicities, and abundance ratios 

SSP-equivalent parameters are derived for each G93 galaxy 
by choosing, for each model 1^, the best-fitting age f , metal- 
licity [Z/H], and enhancement ratio [E/Fe]. Solving for three 

"^We have ignored the dependence of the Hne strength indices on Ti, as TB95 make contradictory statements about its inclusion in their model atmospheres. 
Although their tables include the effects of varying Ti, they clearly state that they have not included TiO lines in their line lists. This will affect the line strengths in 
the coolest giants. However, H/3. Mgb, and (Fe) are little affected by Ti in their models; see their Tables 4-6. 

^This equation assumes that the percentage index change is constant for each step of 0.3 dex in abundance. This assures that index values approach zero gracefully 
at low abundances but predicts infinite indices at high abundances, which is impossible. The scaling law should therefore probably not be applied at levels much 
above [Xi/H] =+0.6. 



Tables ^ and ^ give derived (f, [Z/H], [E/Fe]) values 
and associated uncertainties in the r^/S and rg/2 aper- 
tures, respectively. Errors were derived by searching the 
grid at iHp±<jH(3,Mgb,{Fe)), (H/3,Mg/7 ± ^Mgi,, (Fe)), and 
(H(3,Mgb,{Fe) ± cr^pe)) and taking the maximum deviations 
max(Af), max(A[Z/H]), and max(A[E/Fe]) as the associated 
uncertainties.'' 

The derived SSP-equivalent parameters should be treated 
with caution for the extrapolated solutions (f > 18 Gyr at all 
metaflicities, [Z/H] > 0.5 at aU ages, and t < 8 Gyr at [Z/H] < 
-0.225). However, in the re/8 aperture, which we concentrate 
on in this and the following paper, only one galaxy (NGC 5813) 
has f > 18 Gyr, and only a few more have [Z/H] > 0.5 for any 
enhancement model. The extrapolations are more significant 
for the stellar population parameters in the r^/l aperture (Ta- 
ble ^b|). However, many of these would also lessen or disap- 
pear if TWFBG98 velocity corrections to H/3 were substituted 
for those of G93, or if small H/3 emission fill-in errors were 
corrected (see Appendix^. 

We have checked our fitting procedure by correcting the ob- 
served line strengths back to solar abundance ratios using the 
TB95 response functions for the solved-for values of [Z /H] and 
[E/Fe]. The resulting corrected line strengths are presented in 
Figures I (for r,/8) and g (for re/2) with the W94 models over- 
plotted. These are the predicted line strengths that would be 
seen if the populations had the same f and [Z/H] but [E/Fe] 
= 0. Metallicities and ages inferred from Mgb and (Fe) now 
agree, suggesting that our method for finding for the best-fitting 
parameters by searching in the three-dimensional grid is work- 
ing correctly. These corrected points show graphically our final 
values of t and [Z/H]. 



These errors faithfully . 



agnitude of the uncertainties but not their correlations. Correlated errors in [Z/H] and t can be important, driven jointly by 



observational errors in H/3 (ITrager 199'7|). Fortunately the G93 errors are so small that observationally driven correlations in the output parameters are not important. 
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Derived stellar parameters from the four enhancement mod- 
els are compared in Figure ||. The most notable difference is be- 
tween model 3 (C down, O up) versus all other models: galaxies 
are older, more metal-poor, and less enhanced in model 3 than 
in the others. These differences are driven entirely by the low 
C abundance in model 3; reducing C increases Mgb but has lit- 
tle effect on (Fe). Models with low C (like model 3) therefore 
result in lower overall metallicities, smaller [E/Fe], and older 
ages, as may be seen by following through the consequences of 
a higher Mg b response function in Figure |] 

Is model 3 in fact compatible with observed galaxy line 
strengths? To test this, we augment the G93 indices with data 
on the C-sensitive feature C24668 from the Lick/IDS sample of 
TWFBG98. For each population model, we use the response 
functions of TB95 to compute predicted line strengths for three 
new features — C24668, Mgi, and Mg2 — none of which were 
used in the original fits. Observed versus predicted indices are 
shown in Figure ^. Enhancement model 3, in which C is de- 
pressed, clearly fails systematically to reproduce the strengths 
of the new indices, especially C24668. Models 1, 2, and 4 
are nearly indistinguishable, as expected since the C abundance 
hardly varies among them (cf. Table Model 4 is marginally 
the best (~ la) on account of its slightly higher C abundance, 
a further slight boost for our preferred model. Although model 
4 fits best, it still fails systematically to reproduce the highest 
values of C24668, Mgi and, especially, Mg2. This may indicate 
that C (and perhaps Mg) are actually over-enhanced compared 
to the E-group generally and may signal a breakdown in our as- 
sumption that all E-group elements scale in lockstep. Specific 
element abundance ratios will be explored using the full set of 
Lick indices in future papers. 
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Table 6a 

Ages, metallicities and enhancement ratios through the central r(,/8 aperture 



Name 


Model 


Age (Gyr) 


[Z/H] 


[E/Fe] 


[Fe/H] 


[E/H] 


NGC 221 


1 


3.0±0.7 


0.01 ±0.04 


-0.08 ±0.01 


0.08 ± 0.04 


0.00 ± 0.04 




2 


3.0 ±0.5 


0.02 ±0.04 


-0.08 ±0.01 


0.08 ± 0.04 


0.00 ± 0.04 




3 


2.9 ±0.5 


0.07 ±0.04 


-0.06 ±0.01 


0.12 ±0.04 


0.06 ±0.04 




4 


3.0 ±0.7 


0.00 ±0.05 


-0.08 ±0.01 


0.07 ±0.05 


-0.01 ±0.05 


NGC224 


1 


6.2 ±1.5 


0.37 ±0.06 


0.18±0.02 


0.21 ±0.06 


0.39 ±0.06 




2 


6.4 ± 1.6 


0.35 ±0.05 


0.18 ±0.02 


0.21 ±0.05 


0.39 ±0.05 




3 


6.7 ±1.5 


0.29 ±0.04 


0.12±0.01 


0.20 ±0.04 


0.32 ±0.04 




4 


6.1 ±1.5 


0.38 ±0.07 


0.19 ±0.02 


0.20 ±0.07 


0.39 ±0.07 


NGC 315 


1 


5.8±1.1 


0.31 ±0.04 


0.23 ±0.02 


0.10 ±0.04 


0.33 ±0.04 




2 


6.1 ±1.3 


0.28 ±0.05 


0.23 ±0.02 


0.10 ±0.05 


0.33 ±0.05 




3 


6.7 ±1.5 


0.21 ±0.04 


0.16 ±0.01 


0.09 ±0.04 


0.25 ±0.04 




4 


5.8 ±1.2 


0.32 ±0.06 


0.24 ±0.02 


0.10 ±0.06 


0.34 ±0.06 


NGC 507 


1 


7.6 ±2.9 


0.17 ±0.07 


0.19 ±0.03 


0.00 ±0.08 


0.19 ±0.07 




2 


7.8 ±3.0 


0.15 ±0.07 


0.19 ±0.02 


0.00 ±0.07 


0.19 ±0.07 




3 


9.0 ±2.3 


0.07 ±0.07 


0.13 ±0.02 


-0.03 ±0.07 


0.10 ±0.07 




4 


7.5 ±2.5 


0.18 ±0.06 


0.20 ±0.03 


-0.01 ±0.07 


0.19 ±0.06 


NGC 547 


1 


9.9 ±2.3 


0.20 ±0.05 


0.25 ±0.01 


-0.03 ±0.05 


0.22 ±0.05 




2 


10.2 ±2.5 


0.18 ±0.05 


0.25 ± 0.02 


-0.01 ±0.05 


0.24 ± 0.05 




3 


11.6 ±2.7 


0.09 ±0.05 


0.17 ±0.01 


-0.04 ±0.05 


0.13 ±0.05 




4 


9.5 ±2.2 


0.22 ±0.05 


0.26 ±0.01 


-0.02 ±0.05 


0.24 ±0.05 


NGC 584 


1 


2.6 ±0.3 


0.46 ±0.03 


0.21 ±0.01 


0.27 ±0.03 


0.48 ±0.03 




2 


2.7 ±0.3 


0.43 ±0.03 


0.21 ±0.01 


0.27 ±0.03 


0.48 ±0.03 




3 


2.8 ±0.4 


0.36 ±0.03 


0.14 ±0.01 


0.25 ±0.03 


0.39 ±0.03 




4 


2.5 ±0.3 


0.48 ±0.03 


0.22 ±0.01 


0.28 ±0.03 


0.50 ±0.03 


NGC 636 


1 


4.0 ±0.4 


0.35 ±0.04 


0.11 ±0.01 


0.25 ± 0.04 


0.36 ±0.04 




2 


4.0±0.6 


0.34 ±0.05 


0.11 ±0.02 


0.25 ±0.05 


0.36 ±0.05 




3 


4.3 ±0.5 


0.28 ±0.04 


0.07 ±0.01 


0.23 ±0.04 


0.30 ±0.04 




4 


4.1 ±0.7 


0.34 ±0.07 


0.11 ±0.02 


0.24 ±0.07 


0.35 ±0.07 


NGC 720 


1 


4.7 ±2.4 


0.41 ±0.14 


0.31 ±0.04 


0.13±0.14 


0.44 ±0.14 




2 


4.9 ±2.5 


0.36±0.14 


0.31 ±0.04 


0.12±0.14 


0.43 ±0.14 




3 


6.0±2.4 


0.25 ±0.10 


0.20 ±0.03 


0.10 ±0.10 


0.30 ±0.10 




4 


4.5 ±2.7 


0.44 ±0.15 


0.33 ±0.04 


0.13 ±0.15 


0.46 ±0.15 


NGC 821 


1 


7.8± 1.5 


0.21 ±0.03 


0.14 ±0.01 


0.08 ± 0.03 


0.22 ±0.03 




2 


8.2± 1.0 


0.19±0.03 


0.14 ±0.01 


0.08 ± 0.03 


0.22 ±0.03 




3 


9.0±1.3 


0.14±0.03 


0.10 ±0.01 


0.06 ± 0.03 


0.16 ±0.03 




4 


7.7±1.3 


0.22 ±0.03 


0.15 ±0.01 


0.08 ±0.03 


0.23 ±0.03 


NGC 1453 


1 


7.9 ± 1.4 


0.29 ±0.04 


0.20 ±0.02 


0. 1 1 ± 0.04 


0.31 ±0.04 




2 


8.2± 1.8 


0.26 ± 0.04 


0.20 ±0.02 


0. 1 1 ± 0.04 


0.31 ±0.04 




3 


9.4 ±2.0 


0.19 ±0.04 


0.14 ±0.01 


0.08 ± 0.04 


0.22 ±0.04 




4 


7.9±1.4 


0.30 ±0.04 


0.21 ±0.02 


0.10 ±0.04 


0.31 ±0.04 


NGC 1600 


1 


8.5 ±1.7 


0.34 ± 0.05 


0.21 ±0.02 


0.15 ±0.05 


0.36 ±0.05 




2 


8.9±1.7 


0.31 ±0.05 


0.21 ±0.02 


0.15 ±0.05 


0.36 ±0.05 




3 


9.3 ±2.2 


0.24 ±0.05 


0.14 ±0.01 


0.13 ±0.05 


0.27 ±0.05 




4 


8.6±1.7 


0.35 ±0.05 


0.22 ±0.02 


0.15 ±0.05 


0.37 ±0.05 


NGC 1700 


1 


2.3 ±0.3 


0.49 ±0.03 


0.16 ±0.01 


0.34 ±0.03 


0.50 ±0.03 




2 


2.4 ±0.4 


0.46 ±0.04 


0.16 ±0.01 


0.34 ±0.04 


0.50 ±0.04 




3 


2.6 ±0.4 


0.40 ± 0.04 


0.10 ±0.01 


0.32 ±0.04 


0.42 ±0.04 




4 


2.3 ±0.3 


0.50 ± 0.03 


0.16 ±0.01 


0.35 ± 0.03 


0.51 ±0.03 


NGC 2300 


1 


6.5 ±1.6 


0.34 ±0.05 


0.23 ±0.02 


0.13 ±0.05 


0.36 ±0.05 




2 


6.8±1.8 


0.31 ±0.05 


0.23 ±0.02 


0.13 ±0.05 


0.36 ±0.05 




3 


7.2±1.5 


0.24 ±0.04 


0.15±0.01 


0.13 ±0.04 


0.28 ±0.04 




4 


6.3 ±1.4 


0.36 ±0.04 


0.24 ±0.02 


0.14 ±0.04 


0.38 ±0.04 


NGC 2778 


1 


5.4 ±2.1 


0.28 ±0.06 


0.22 ±0.02 


0.08 ±0.06 


0.30 ±0.06 




2 


5.8±1.9 


0.25 ±0.06 


0.22 ±0.02 


0.08 ±0.06 


0.30 ±0.06 




3 


6.8±1.8 


0.18 ±0.05 


0.15 ±0.01 


0.07 ±0.05 


0.22 ±0.05 




4 


5.4±1.7 


0.29 ±0.07 


0.23 ±0.02 


0.08 ±0.07 


0.31 ±0.07 



TRACER ETAL. 
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Name 


Model 


Age (Gyr) 


[Z/H] 


[E/Fe] 


[Fe/H] 


[E/H] 


NGC 3377 


1 


3.7 ±0.9 


0.19 ±0.05 


A 1 A 1 A A'^ 

0.19±0.02 


A A'^ 1 A Ag 

0.02 ± 0.05 


A A 1 1 A Ag 

0.21 ±0.05 




o 

L 


4.1 it U.o 


U.i J it U.Uj 


u.iy it u.ui 


A AA -U A A'^ 
U.UU it U.Uj 


A 1 Q -U A (\'\ 
U.iy it U.Uj 




■3 


4.D it U.J 


n no -J- n n7 
u.uy it u.uz 


u.i J it U.Ul 


_A AI -U A A7 
—U.Ul it U.UZ 


A 1 7 -U A A7 
U.IZ it U.UZ 




/I 

4 


j.o it u.y 


u.iy it u.uo 


f\ 1 Q _i_ f\ n7 
U.iy it u.uz 


A AI J_ A A/; 
U.Ul it U.UO 


A 7A _!_ A A/; 
U.ZU it U.UO 


NGC 3379 


1 


8.9di 1.9 


U.zU ± 0.04 


A ^ A f A A 1 

0.20 ± 0.01 


A A'^ 1 A A ^ 

0.02 ± 0.04 


A 1 A A /I 

0.22 ±0.04 




z 


O -7 1 1 O 

y. / it i.y 


U. i / it U.U4 


n 7n -1- n ni 
U.ZU it U.Ul 


A n7 -1- A A/l 

U.UZ zt U.U4 


A 77 -I- A A/1 
U.ZZ it U.U4 






iU.O it i.D 


U. iU zt U.Uj 


U. 14 it U.Ul 


A n 1 -1- A ni 
—U.U 1 zt U.Uj 


A 1 Q -1- A f^l 
U. 1 J zt U.UJ 






Q c _|_ 1 
o.S it i.D 


U . Z i It U . U4 


n 7 1 -4- f 1 n 1 

U.Zi It VJ.U 1 


U.U i zt U.U4 


A 77 -1- A (\A 
U.ZZ zt U.U4 


NGC 3608 


1 


6.9 di 1.5 


0.25 ±0.04 


0.16 ±0.02 


0.10 ±0.04 


0.26 ±0.04 




L 


1 % -X- \ 
1 .j± i.O 


U.Zj It U.U4 


n 1 -1- n n7 
u. 10 It u.uz 


A 1 1 -I- A A/1 
U. 1 1 It U.U4 


A 77 -I- A A/1 
U.Z / It U.U4 




•3 


Q n -I- 1 Q 

o.U it i.O 


^"1 1 7 _l_ ^"1 f\A 
u.i/ It U . U4 


U . 1 i It u . U i 


A AQ -I- A A/1 
U.uy It U.U4 


A 7A -I- A A/1 
U.ZU It U.U4 




4 


O.O it i.J 


U.ZO lb U.U4 


n 1 7 j_ n n7 
U. 1 / it U.UZ 


A 1 A _!_ A A/1 
U. lU it U.U4 


A 77 _!_ A A/1 
U.Z / it U.U4 


NGC 3818 


1 


5.9 lb 1.5 


0.34 ± 0.07 


0.22 ±0.02 


0.14 ±0.07 


0.36 ±0.07 




L 


o.U it i.j 


U. jZ zt U.Uj 


n 77 _L n m 
U.ZZ it U.UZ 


A 1 ^ -1- A A^ 

U. 1 J zt U.Uj 


A 17 _L A A^ 

U.J / it U.Uj 




'2 


O.O it i . / 


n 7/1 J- n 
U.Z4 it vJ.Uj 


U. 1 J it U.Ul 


A 1 1 _!_ A A^ 

U. 1 J it U.Uj 


A 7Q _!_ A A^ 

U.Zo it U.Uj 




4 


J .0 It i.O 


u. jO It u.uo 


n 71 _i_ n m 
U.Zj it U.UZ 


U. i J it U.UO 


A IS -I- A A*^ 
U.Jo it u.uo 


NGC 4261 


1 


16.1 ±3.2 


0.17 ±0.04 


0.19 ±0.01 


A AA 1 A A ^ 

0.00 ±0.04 


A -1 A 1 f\ r\ A 

0.19 ±0.04 




L 


iO.O it Z.J 


U. i J it U.Uj 


U. ly it U.Ul 


A AA -I- A AI 
U.UU it U.Uj 


A 1 Q _!_ A AI 

u. ly it U.UJ 






i / .y it z.o 


u.uy It U.Uj 


U . 1 J It u , u i 


A A 1 -I- A AI 
-U.Ul It U.UJ 


A 1 7 _l_ A AI 
U. iZ It U.UJ 




A 

4 


ij.o it z.y 


U. io it U.U4 


n 7n -U /^^ n 1 
U.ZvJ zt U.Ul 


A A 1 _!_ A A/1 
—U.Ul it U.U4 


A 1 Q -U A A/1 

u. ly it U.U4 


NGC 4374 


1 


13. 2± 1.6 


0.10 ±0.03 


0.20 ±0.01 


-0.08 ± 0.03 


0.12 ±0.03 




L 


Lj./ It i.O 


n HQ -I- n m 

U.Uo It U.Uj 


U.ZVJ zt VJ.Ui 


A A7 -I- A AI 
—U.U / It U.UJ 


A 1 1 _l_ A AI 
U. 1 J it U.UJ 




'1 
D 


1 < 1 _L 1 /; 
i D. i it i .0 


U.UU it U.Uj 


U. 14 zt U.Ul 


A 1 1 _!_ A AI 
— U. 1 1 it U.UJ 


A AI _!_ A AI 
U.UJ it U.UJ 




4 


iZ. / it z.u 


U. iZ it U.Uj 


n 7n -U n ni 
U.ZU it U.Ul 


A A7 _!_ A AI 
—U.U / it U.UJ 


A 1 1 _!_ A AI 
U. 1 J it U.UJ 


NGC 4472 


1 


8.0 ±2.0 


0.24 ± 0.04 


0.20 ± 0.02 


0.06 ± 0.04 


0.26 ±0.04 




L 


O.J it Z.i 


U.ZZ it U.U4 


n 7 1 J- n n7 
U.ZI it U.UZ 


A A/^ _!_ A A/1 
U.UO it U.U4 


A 77 _!_ A A/1 
U.Z / it U.U4 




-2 
J 


1 n 1 -U 

iU. i It z.z 


f\ ^ A -\- r\ f\A 
U. 14 it U.U4 


U. 14 it U.Ul 


A AI -I- A A/1 
U.UJ it U.U4 


A 1 7 _l_ A A/1 
U. i / It U.U4 




4 


/.y it z.i 


n 7'\ -U n f\< 
U.Zj it U.Uj 


n 7 1 -U n n7 
U.ZI it U.uz 


A A^ -U A A^ 
U.Uj it U.Uj 


A 7/^ -U A A^ 
U.ZO it U.Uj 


NGC 4478 


1 


4.7 ±2.0 


0.28 ±0.09 


f\ -t A 1 A A^ 

0.14 ±0.03 


A i g 1 A AA 

0.15 ±0.09 


A AA 1 A AA 

0.29 ±0.09 




L 


4.0 It Z.J 


f\ 7A -1- n no 
U.ZO It u.uy 


U. 14 It U.Uj 


A 1 '^ -I- A AQ 
U. 1 J It U.uy 


A 70 -U A AQ 

U.zy It U.uy 




I 

D 


J.J it Z.i 


U.Zi it U.U/ 


n no -1- (\ m 
U.uy it U.uz 


U. 14 zt U.U / 


A 71 _L- A l^n 

U.ZJ it U.U / 




4 


4.D it Z.J 


n 70 1 n 1 n 

u.zy it u.iu 


U.l J it U.Uj 


A 1 ^ _!_ A 1 A 

U.lj it U.IU 


A 1A 1 A 1 A 
U.jU it U.IU 


NGC 4489 


1 


2.5 ±0.4 


A 1 ^ 1 A A^ 

0.14 ±0.06 


A AI 1 A Ai^ 

0.03 ± 0.02 


A 1 i 1 A A/T 

0.11 ±0.06 


Aid 1 A A/T 

0.14 ±0.06 




L 


Z.J It U.4 


U.14 it U.UO 


n A7 _i_ n n7 
U.UZ it U.UZ 


A 1 7 _1_ A A/i 

U.IZ it U.UO 


A 1 /I J_ A C\t^ 

U.14 it U.UO 




'1 
D 


Z.J It u.4 


U.I J it U.UO 


n A7 _i_ n A7 
U.UZ it U.UZ 


A 1 1 _L- A A/i 

U.ll it U.UO 


A 11 J_ A C\f^ 

U.l J it u.uo 




A 
4 


Z.J It U.4 


U.14 it U.uo 


U.Uj it U.uz 


A 1 1 -1- A Afi 
U.ll dl U.UO 


A 1 4 -U A Aft 
U.14 it U.UO 


NGC 4552 


1 


10.6 ± 1.5 


/\ f^C 1 A A J 

0.25 ±0.04 


A /^/^ 1 A A i 

0.22 ±0.01 


A Ag 1 A A *i 

0.05 ± 0.04 


A A^ 1 A A ji 

0.27 ±0.04 




L 


11 A -U 7 1 
1 1.4 It Z. 1 


U.ZZ It 


U.ZZ It U.Ul 


U.UD It U.U4 


n 07 4- n (\A 

KJ.ZJ It U.U4 




'X 
o 


1 7 S -U 7 1 
IZ.O It Z. 1 


U.I J It U.Uj 


u.l J It U.Ul 


U.U4 It U.UJ 


A 1 Q 4_ A 

u.iy It U.Uj 




4 


1 n C _L 1 A 
lU.J It 1.4 


U.Z/ itU.U4 


U.Zo It U.Ul 


u.uo it U.U4 


A OO 4- A f\A 

u.zy It U.U4 


NGC 4649 


1 


11.9 ±2.8 


0.25 ±0.04 


0.23 ±0.01 


0.04 ±0.04 


0.27 ±0.04 




2 


12.7 ±2.3 


0.23 ±0.03 


0.23 ±0.01 


0.05 ±0.03 


0.28 ±0.03 




3 


15.0 ±2.9 


0.15 ±0.03 


0.16±0.01 


0.03 ±0.03 


0.19 ±0.03 




4 


11.9±2.6 


0.27 ±0.04 


0.24 ±0.01 


0.05 ±0.04 


0.29 ±0.04 


NGC 4697 


1 


8.9 ±2.0 


0.05 ±0.05 


0.10 ±0.01 


-0.04 ±0.05 


0.06 ±0.05 




2 


9.1 ±2.0 


0.04 ±0.05 


0.10 ±0.01 


-0.04 ±0.05 


0.06 ±0.05 




3 


9.5 ±1.9 


0.00 ±0.05 


0.07 ±0.01 


-0.05 ±0.05 


0.02 ±0.05 




4 


8.8±2.0 


0.06 ±0.05 


0.10 ±0.02 


-0.03 ±0.05 


0.07 ± 0.05 


NGC 5638 


1 


8.7±1.3 


0.18±0.03 


0.18 ±0.01 


0.02 ± 0.03 


0.20 ±0.03 




2 


9.5±1.5 


0.15 ±0.03 


0.18±0.01 


0.01 ±0.03 


0.19 ±0.03 




3 


10.1 ±1.4 


0.09 ±0.03 


0.13±0.01 


-0.01 ±0.03 


0.12 ±0.03 




4 


8.5±1.5 


0.19 ±0.03 


0.19±0.01 


0.01 ±0.03 


0.20 ± 0.03 


NGC 5812 


1 


5.5±1.0 


0.37 ±0.03 


0.19 ±0.01 


0.20 ±0.03 


0.39 ± 0.03 




2 


5.9±1.0 


0.34 ±0.03 


0.19±0.01 


0.19 ±0.03 


0.38 ±0.03 




3 


6.7±1.1 


0.27 ±0.03 


0.12 ±0.01 


0.18±0.03 


0.30 ±0.03 




4 


5.4 ±1.0 


0.38 ±0.04 


0.20 ±0.01 


0.19 ±0.04 


0.39 ±0.04 



STELLAR POPULATIONS OF EARLY-TYPE GALAXIES. I. 
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Name 


Model 


Age (Gyr) 


[Z/H] 


[E/Fe] 


[Fe/H] 


[E/H] 




\ 




W.UJ HI L/.U9 


L/.ZL/ ni L/.Ul 


VJ.Z9 It \J.\JJ 


V.KJJ It \J.\JJ 




2 


18.9 ±2.4 


-0.07 ± 0.03 


0.21 ±0.02 


-0.23 ± 0.03 


-0.02 ± 0.03 




3 


19.4 ±2.0 


-0.12 ±0.03 


0.14 ±0.01 


-0.23 ± 0.03 


-0.09 ± 0.03 




4 


18.5 ±2.3 


-0.04 ±0.03 


0.21 ±0.01 


-0.24 ± 0.03 


-0.03 ± 0.03 


iNvJ\w- Do J 1 


1 
1 




u . 9 i znyi.yfj 


U. io It V1.V7 1 


-A- OS 
\f.JJ ZIZ V7.V19 


S^ ^ OS 
VI. 9 9 ZIZ V1.V79 




2 


2.7 ±0.3 


0.50 ±0.03 


0.18±0.01 


0.36 ±0.03 


0.54 ±0.03 




3 


2.9 ±0.5 


0.43 ± 0.05 


0.12 ±0.01 


0.34 ±0.05 


0.46 ± 0.05 




4 


2.6 ±0.3 


0.54 ±0.04 


0.19±0.01 


0.36 ±0.04 


0.55 ± 0.04 


iNvJV. JofD 


1 
i 




U. iZ ^ U.U9 


n 90 -1- 09 
U.ZU It U.UZ 


0(^ -1- OS 
— V1.VJO It VJ.VJ9 


1 4 4- OS 
VJ. i't It U.U9 




2 


14.9±3.3 


10 + 05 


0.21 ±0.01 


-0.06 ± 0.05 


1 5 ± 05 




3 


16.2 ±2.9 


0.03 ±0.04 


0.14 ±0.01 


-0.08 ± 0.04 


0.06 ± 0.04 




4 


14.1 ±2.9 


0.14 ±0.05 


0.21 ±0.02 


-0.06 ± 0.05 


0.15 ±0.05 




^ 


11 Q -1- 9 Q 


u.iuniu.u^ 


99 4-0 01 

U.ZZ ^ U.U i 


—0 (^A -t- (^A 


18 4-0 04 

VI . i O 31 VJ . 




2 


12.9 ±2.0 


0.13 ±0.03 


0.22 ±0.02 


-0.04 ± 0.03 


0.18 ±0.03 




3 


14.7 ±2.6 


0.05 ±0.04 


0.15±0.01 


-0.06 ± 0.04 


0.09 ±0.04 




4 


11.8±2.3 


0.17 ±0.03 


0.23 ±0.01 


-0.04 ± 0.03 


0.19 ±0.03 






1 s -t- n 9 


\j.uy zn \). 1 1 


1 s + f1 0^ 
u. 1 J zn VI. u J 


55 ± 11 


VI. /VI ZIZ VI. 1 1 




2 


1.5 ±0.2 


0.68 ±0.11 


0.15 ±0.02 


0.56 ±0.11 


0.71 ±0.11 




3 


1.6 ±0.1 


0.58 ±0.05 


0.10 ±0.01 


0.50 ±0.05 


0.60 ± 0.05 




4 


1.5 ±0.1 


0.70 ± 0.06 


0.15 ±0.03 


0.56 ± 0.07 


0.71 ±0.06 




1 
1 


^.j m u.o 


U. 911 ZJZ VI. UD 


1 S ^ 09 
VI. 19 ZIZ VI. VIZ 


1^-1-0 Oft 
VI. iO It V1.V10 


0^1 .^0 Oft 
VI. 9 i ZIZ V1.V1D 




2 


4.4 ±0.9 


0.28 ± 0.05 


0.15 ±0.02 


0.16 ±0.05 


0.31 ±0.05 




3 


4.8± 1.2 


0.22 ±0.06 


0.10 ±0.01 


0.14 ±0.06 


0.24 ± 0.06 




4 


4.4 ±0.7 


0.30 ±0.06 


0.15±0.02 


0.16 ±0.06 


0.31 ±0.06 


iNvJV. /UJZ 


1 


1 9 Q -t- 4 




9'^ -t- 01 
U.Z9 It U.V1I 


oft 4- OS 
— V1.VJO It V1.VJ9 


1 7 ^ OS 
VI. i / zn v.yjj 




2 


13.9±2.9 


0.12 ±0.05 


0.23 ±0.02 


-0.06 ± 0.05 


0.17 ±0.05 




3 


15.3±3.2 


04 + 04 


0.16±0.01 


-0.08 ± 0.04 


0.08 ± 0.04 




4 


12.7 ±3.0 


0.16 ±0.05 


0.24 ±0.01 


-0.06 ± 0.05 


0.18 ±0.05 


NGC 7454 


1 


5.1 ±0.9 


-0.07 ± 0.03 


0.06 ±0.02 


-0.12 ±0.04 


-0.06 ± 0.03 




2 


5.1±1.0 


-0.07 ± 0.04 


0.06 ±0.02 


-0.12 ±0.04 


-0.06 ± 0.04 




3 


5.2±1.0 


-0.09 ±0.03 


0.04 ±0.01 


-0.12 ±0.03 


-0.08 ±0.03 




4 


5.0±1.0 


-0.06 ±0.04 


0.06 ±0.02 


-0.12 ±0.04 


-0.06 ±0.04 


NGC 7562 


1 


7.9 ±1.3 


0.19 ±0.03 


0.16 ±0.01 


0.04 ±0.03 


0.20 ±0.03 




2 


8.2±1.4 


0.17 ±0.04 


0.16 ±0.02 


0.05 ±0.04 


0.21 ±0.04 




3 


9.3 ±1.7 


0.11 ±0.04 


0.11 ±0.01 


0.03 ±0.04 


0.14 ±0.04 




4 


7.7 ±1.4 


0.20 ±0.04 


0.17 ±0.01 


0.04 ±0.04 


0.21 ±0.04 


NGC 7619 


1 


15.1 ±2.6 


0.19 ±0.03 


0.17 ±0.01 


0.03 ±0.03 


0.20 ±0.03 




2 


15.8 ±2.4 


0.17 ±0.03 


0.17 ±0.01 


0.04 ±0.03 


0.21 ±0.03 




3 


16.5 ±1.5 


0.12 ±0.03 


0.12 ±0.01 


0.03 ±0.03 


0.15 ±0.03 




4 


14.8 ±2.3 


0.20 ±0.03 


0.18 ±0.01 


0.03 ±0.03 


0.21 ±0.03 


NGC 7626 


1 


13.6 ±2.3 


0.14 ±0.03 


0.24 ±0.01 


-0.08 ±0.03 


0.16 ±0.03 




2 


14.6 ±2.6 


0.11 ±0.04 


0.24 ±0.01 


-0.08 ±0.04 


0.16 ±0.04 




3 


16.2 ±2.4 


0.03 ±0.03 


0.16 ±0.01 


-0.09 ±0.03 


0.07 ±0.03 




4 


12.9 ±2.9 


0.16 ±0.04 


0.25 ±0.01 


-0.07 ±0.04 


0.18 ±0.04 


NGC 7785 


1 


8.8 ±1.9 


0.19 ±0.04 


0.16 ±0.01 


0.04 ±0.04 


0.20 ±0.04 




2 


9.2 ±1.9 


0.17 ±0.04 


0.16 ±0.01 


0.05 ±0.04 


0.21 ±0.04 




3 


9.9 ±1.9 


0.12 ±0.04 


0.11 ±0.01 


0.04 ±0.04 


0.15 ±0.04 




4 


8.7 ±1.9 


0.20 ±0.04 


0.16 ±0.01 


0.05 ±0.04 


0.21 ±0.04 



TRACER ETAL. 



Table 6b 

Ages, metallicities and enhancement ratios through the global re/2 aperture 



Name 


Model 


Age (Gyr) 


[Z/H] 


[E/Fe] 


[Fe/H] 


[E/H] 


NGC 221 


1 


4.9± 1.1 


-0.08 ± 0.04 


-0.07 ± 0.02 


-0.02 ± 0.04 


-0.09 ± 0.04 




2 


4.9± 1.0 


-0.07 ± 0.04 


-0.07 ±0.02 


-0.02 ± 0.04 


-0.09 ± 0.04 




3 


4.9 ±1.0 


-0.06 ±0.03 


-0.05 ±0.01 


-0.02 ±0.03 


-0.07 ±0.03 




4 


4.9±1.3 


-0.08 ±0.05 


-0.07 ±0.02 


-0.01 ±0.05 


-0.08 ±0.05 


NGC 315 


1 


6.9 ±2.4 


0.16±0.07 


0.23 ±0.02 


-0.05 ±0.07 


0.18±0.07 




2 


7.3 ±2.2 


0.13±0.06 


0.23 ±0.02 


-0.05 ± 0.06 


0.18 ±0.06 




3 


8.3 ±2.4 


0.04 ±0.06 


0.15 ±0.02 


-0.07 ±0.06 


0.08 ±0.06 




4 


6.9 ±1.9 


0.17 ±0.06 


0.24 ±0.02 


-0.05 ±0.06 


0.19 ±0.06 


NGC 507 


1 


3.5 ±3.2 


0.27 ±0.15 


0.37 ±0.04 


-0.07 ±0.15 


0.30 ±0.15 




2 


4.1 ±3.1 


0.20 ±0.13 


0.37 ±0.04 


-0.09 ±0.13 


0.28 ±0.13 




3 


4.9 ±3.0 


0.08 ±0.12 


0.25 ±0.02 


-0.11 ±0.12 


0.14 ±0.12 




4 


3.5 ±2.7 


0.28 ±0.13 


0.39 ±0.04 


-0.08 ±0.14 


0.31 ±0.13 


NGC 547 


1 


18.7 ±3.8 


-0.04 ±0.06 


0.24 ±0.02 


-0.26 ±0.06 


-0.02 ±0.06 




2 


18.6±3.3 


-0.05 ±0.04 


0.24 ±0.02 


-0.24 ±0.04 


0.00 ±0.04 




3 


19.7 ±3.2 


-0.13 ±0.05 


0.17 ±0.02 


-0.26 ±0.05 


-0.09 ±0.05 




4 


18.1 ±3.2 


-0.02 ±0.05 


0.25 ±0.02 


-0.25 ±0.05 


0.00 ±0.05 


NGC 584 


1 


3.4 ± 1.2 


0.25 ± 0.07 


0.20 ± 0.02 


0.07 ± 0.07 


0.27 ± 0.07 




2 


3.8 ±0.9 


0.21 ±0.06 


0.20 ± 0.02 


0.06 ± 0.06 


0.26 ± 0.06 




3 


4.4±1.1 


0.13 ±0.06 


0.14 ±0.02 


0.02 ±0.06 


0.16 ±0.06 




4 


3.4±1.1 


0.26 ±0.07 


0.21 ±0.02 


0.06 ±0.07 


0.27 ±0.07 


NGC 636 


1 


6.7 ±1.3 


0.10 ±0.06 


0.11 ±0.02 


0.00 ±0.06 


0.11 ±0.06 




2 


6.8 ±1.3 


0.09 ±0.06 


0.11 ±0.02 


0.00 ±0.06 


0.11 ±0.06 




3 


7.2±1.3 


0.05 ±0.05 


0.07 ±0.01 


0.00 ±0.05 


0.07 ±0.05 




4 


6.8±1.4 


0.10 ±0.05 


0.11 ±0.02 


0.00 ±0.05 


0.11 ±0.05 


NGC 720 


1 


1.5 ±0.4 


0.90 ± 0.23 


0.41 ±0.16 


0.53 ±0.27 


0.94 ± 0.23 




2 


1.7 ±0.4 


0.75 ±0.21 


0.40 ±0.08 


0.44 ±0.22 


0.84±0.21 




3 


1.8 ±0.7 


0.59 ±0.15 


0.27 ±0.03 


0.39±0.15 


0.66 ±0.15 




4 


1.1 ±0.8 


1.13 ±0.42 


0.56±0.18 


0.61 ±0.45 


1.17±0.42 


NGC 821 


1 


7.4± 1.4 


0.10±0.05 


0.12 ±0.02 


-0.01 ±0.05 


0.11 ±0.05 




2 


7.6± 1.4 


0.08 ±0.05 


0.12 ±0.02 


-0.01 ±0.05 


0.11 ±0.05 




3 


7.9 ±1.4 


0.04 ± 0.04 


0.08 ± 0.02 


-0.02 ± 0.04 


0.06 ±0.04 




4 


7.1 ±1.2 


0.12 ±0.05 


0.13 ±0.02 


0.00 ±0.05 


0.13 ±0.05 


NGC 1453 


1 


8.7 ±2.6 


0.14 ±0.06 


0.15 ±0.02 


0.00 ± 0.06 


0.15 ±0.06 




2 


9.1±2.5 


0.12 ±0.06 


0.15 ±0.02 


0.00 ± 0.06 


0.15 ±0.06 




3 


9.9 ±2.6 


0.07 ±0.06 


0.10 ±0.02 


-0.01 ±0.06 


0.09 ± 0.06 




4 


8.5 ±2.7 


0.15 ±0.07 


0.16 ±0.02 


0.00 ± 0.07 


0.16 ±0.07 


NGC 1600 


1 


4.4± 1.8 


0.50 ±0.14 


0.28 ±0.03 


0.24 ±0.14 


0.52 ±0.14 




2 


4.6± 1.9 


0.46 ±0.13 


0.28 ± 0.03 


0.24 ±0.1 3 


0.52±0.13 




3 


5.2± 1.8 


0.35 ±0.08 


0.18 ±0.02 


0.21 ±0.08 


0.39 ±0.08 




4 


4.3 ±1.5 


0.53 ±0.13 


0.29 ±0.04 


0.26 ±0.14 


0.55 ±0.13 


NGC 1700 


1 


2.8 ±0.6 


0.32 ±0.06 


0.13 ±0.02 


0.20 ± 0.06 


0.33 ±0.06 




2 


2.9 ±0.5 


0.29 ±0.04 


0.13 ±0.02 


0.19 ±0.04 


0.32 ±0.04 




3 


3.0±0.9 


0.25 ±0.05 


0.08 ± 0.02 


0.19 ±0.05 


0.27 ±0.05 




4 


2.8 ±0.5 


0.32 ±0.05 


0.13 ±0.02 


0.20 ±0.05 


0.33 ±0.05 


NGC 2300 


1 


10.2±1.9 


0.13 ±0.04 


0.21 ±0.02 


-0.06 ± 0.04 


0.15 ±0.04 




2 


10.6±1.8 


0.10 ±0.04 


0.22 ±0.02 


-0.07 ± 0.04 


0.15±0.04 




3 


11.5±2.3 


0.02 ±0.04 


0.15 ±0.02 


-0.09 ± 0.04 


0.06 ±0.04 




4 


10.1 ±2.0 


0.14 ±0.04 


0.22 ± 0.02 


-0.06 ± 0.04 


0.16 ±0.04 


NGC 2778 


1 


15.1±3.6 


-0.05 ± 0.06 


0.18 ±0.02 


-0.21 ±0.06 


-0.03 ± 0.06 




2 


15.1±3.4 


-0.06 ± 0.05 


0.18 ±0.02 


-0.20 ± 0.05 


-0.02 ± 0.05 




3 


15.7±3.6 


-0.12 ±0.06 


0.13 ±0.02 


-0.22 ±0.06 


-0.09 ±0.06 




4 


14.9 ±3.5 


-0.04 ±0.06 


0.19 ±0.02 


-0.22 ±0.06 


-0.03 ±0.06 


NGC 3377 


1 


6.0 ±1.2 


-0.13 ±0.05 


0.16 ±0.02 


-0.28 ±0.05 


-0.12 ±0.05 




2 


6.1±1.1 


-0.14 ±0.04 


0.16 ±0.02 


-0.26 ±0.04 


-0.10 ±0.04 




3 


6.3±1.3 


-0.18 ±0.04 


0.11 ±0.01 


-0.26 ±0.04 


-0.15 ±0.04 




4 


5.9±1.2 


-0.12 ±0.04 


0.17 ±0.02 


-0.28 ±0.04 


-0.11 ±0.04 



STELLAR POPULATIONS OF EARLY-TYPE GALAXIES. I. 



Table 6b — Continued 



Name 


Model 


Age (Gyr) 


[Z/H] 


[E/Fe] 


[Fe/H] 


[E/H] 


NGC 3379 


1 


13.4 ±2.4 


-0.01 ±0.04 


0.17 ±0.02 


-0.17 ±0.04 


0.00 ± 0.04 




L 


117-1-0 1 
1 J. / it Z.l 


— U.Uj it U.Uj 


U.lcS it U.UZ 


_A 1 7 _|_ A CVX 
—U.l / it U.Uj 


A A1 _|_ A (XX 

U.UI it U.Uj 




'X 
J 


14.J it Z.J 


— U.Uo it U.Uj 


A 1 -U n m 

U.IZ it U.Ui 


_A 1 7 _|_ A (YX. 
—U.l / it U.Uj 


_n n^ -U n (yx 

—U.Uj it U.Uj 




A 

4 


1 ■J O J- /I 
1 j.Z it Z.4 


(\ f\(\ J- (\ {\A 
U.UU it U.U4 


A 1 Q J- A AO 

U.IIS it U.UZ 


A 1 7 J- A A/1 

— U. 1 / it U.U4 


A A1 J- A A/1 

U.UI it U.U4 


NGC 3608 


1 


8.9di2.8 


0.07 ±0.07 


0.07 ± 0.02 


0.01 ±0.07 


0.08 ±0.07 




L 


O 1 IOC 

y.i it Z.J 


U.UO it U.U/ 


n A"? _L n m 
U.U / it U.UZ 


U.UI ± U.U / 


A AQ -L. A A7 
U.Uo zt U.U/ 






y.z it Z.J 


n ri/i J- n n*^ 
U.U4 it U.UO 


n A^ _L n AO 
U.Uj it U.UZ 


(\ (\(\ _|_ (\ AA 
U.UU zt U.Uo 


A A^ -U A AA 

U.Uj zt U.UO 




/I 
4 


y.u it Z.J 


U.U / It U.U / 


n n7 -I- n no 

U.U / It U.UZ 


\).\)\) zt U.U / 


fl (Vl _1_ (\ (Yl 

vJ.U/ zt vJ.U/ 


NGC 3818 


1 


8.2 ±3.0 


0.04 ± 0.09 


0.16 ±0.03 


-0.11 ±0.09 


0.05 ±0.09 




L 


Q -I- 1 O 
O.J it Z.V 


n m -1- no 
u.uz it u.uy 


n 1 7 4- n Ai 
u. i / it U.Uj 


A 1 1 _|_ A AO 

— u. i 1 It U.uy 


n AA -1- n AO 
U.UO it u.uy 




% 
J 


O T -I- /I 

y.z it z.4 


— VJ.Uj it U.U/ 


n 1 1 -I- A AO 
U. i 1 it U.UZ 


A 1 1 _|_ A A7 
— U. i 1 It U.U / 


n An -1- n A7 

U.UU it U.U/ 




A 

4 


o n 1 1 1 
J5.U it J. i 


n -U n no 
U.Uj it U.uy 


n 1 7 _L n A^ 
U. i / it U.Uj 


A 1 1 _L A AO 

— u. i 1 it U.uy 


A AA -L A AO 

U.UO it U.uy 


NGC 4261 


1 


21.4±2.2 


-0.03 ±0.03 


0.19 ±0.02 


-0.20 ±0.04 


-0.01 ±0.03 




L 


Zi .4 it V.I 


n (\A -L n m 
— U.U4 it U.UZ 


U. iy it U.UI 


A 1 O J- A AO 

— u. iy it U.UZ 


A (\(\ -i- (\ AO 
U.UU zt U.UZ 






Zi .5 it i.D 


n no -L n ni 
—U.uy it U.Uj 


U. i J it U.UI 


A 1 O _L A A^ 

— U. iy it U.Uj 


A AA -1- A AQ 
—U.Uo it U.Uj 




4 


Zi .U it V.I 


— U.Ui it U.Uj 


n 1 Q -I- n AO 
u. iy it U.UZ 


A 1 Q -I- A A/1 
— U. iy It U.U4 


U.UU It U.UJ 


NGC 4374 


1 


14.2 ±2.3 


-0.01 ±0.04 


0.18 ±0.02 


-0.17 ±0.04 


0.01 ±0.04 




Z 


1/1 A -i- O A 
14.4 it Z.4 


n ni J- n n/i 
—U.Uj it U.U4 


n 1 O _L A AO 

u. iy it U.UZ 


A 1 Q -L A A/1 
— U. io it U.U4 


A A 1 -U A A/1 
U.Ui it U.U4 




J 


1 J. J It Z.J 


—U.uy It u,U4 


U. i J It U.VJl 


(\ 1 Q -1- A A/1 
— U. i y zt U.U4 


A AA -I- A A/1 
—U.UO It U.U4 




/I 

4 


14.4 it Z.o 


—U.Ui it U.Uj 


("A 1 O _L ("A (^1 1 

u. iy it U.UI 


A 1 O -J- A 

— u. iy zt u.vjj 


U.UU zt U.Uj 


NGC 4472 


1 


8.6±2.9 


0.17 ±0.07 


0.18 ±0.02 


0.01 ±0.07 


0.19 ±0.07 




Z 


O 1 _|_ O T 

V. 1 it z. / 


U. 14 It U.UO 


n 1 Q -I- n A^ 
u. iy it U.Uj 


n A 1 -I- A f^A 
—U.UI It U.VJO 


{\ 1 _l_ / ^ AA 
VJ, 1 zt VJ.UO 




o 


y.y it z. / 


U.Uo it U.UO 


n 1 ^ _i_ n AO 
U. i J it U.UZ 


A AO _L A AA 
—U.UZ it U.Uo 


A 1 1 _U A AA 
U. 1 i zt U.UO 




4 


0.4 it Z. / 


U. lo it U.UO 


n 1 o -J_ n AO 
U. iy it U.UZ 


A AA _L A AA 
U.UU it U.Uo 


A 1 O -U A AA 
U. ly it U.UO 


NGC 4478 


1 


10.5 ± 1.8 


-0.02 ± 0.04 


0.15 ±0.02 


-0.16 ±0.04 


-0.01 ±0.04 




z 


lU.O it 1.0 


—U.Uj it U.Uj 


n 1 ^ _L A AO 
U. i J it U.UZ 


A 1 ^ _!_ A r\i 
— U. i J it U.Uj 


A AA _L A A1 

U.UU it U.Uj 




a 
J 


11 1 _i_ 1 
1 i . 1 it 1.0 


(\ (Vl _L i^l (\A 

— U . U / It U . U4 


n 1 n -I- A AO 
U. iU it U.UZ 


A 1 ^ -1- A A/1 

— U. i J it U.U4 


A n^ -I- A A/l 
—U.Uj It U.U4 




4 


lU.J it 1. / 


—U.UI it U.U4 


A 1 -U fi no 

U.IO it U.UZ 


—U.IO it U.U4 


n nn -u n (\a 

U.UU it U.U4 


NGC 4489 


1 


4.6 ±0.5 


-0.15 ±0.04 


-0.04 ±0.03 


-0.11 ±0.05 


-0.15±0.04 




o 

Z 


4.0 It U.J 


— U, 1 J It U,U4 


A (\A J~ A (\A 
— U,U4 It U,U4 


A 1 -I- n n'^ 

— U. i Z It U.Uj 


A 1 A _|_ n (\A 
— U. i O It U.U4 




■2 


4.0 it U.O 


— U. 14 it U.Uj 


— U.Uj zt U.UZ 


A 1 -1- A n^; 
— U. iZ zt U.Uj 


A 1 ^ -1- A A^ 

— U. i J it U.Uj 




4 


4.0 it U.J 


A 1 C -1- A f\A 

— U.l J it U.U4 


A A/1 -1- A A/1 

— U.U4 it U.U4 


A 1 1 -I- A A^ 

—U.l 1 it U.Uj 


A 1 ^ -I- A A/1 

—U.l J it U.U4 


NOC 4552 


1 


13.3 ±2.3 


0.09 ± 0.04 


0.21 ±0.02 


-0.10 ±0.04 


r» 1 1 1 r» f\A 

0.11 ±0.04 




Z 


ij.y ± Z.J 


A -J- A A/1 

U.UO ± U.U4 


A OO -L A AO 
U.ZZ it U.UZ 


A 1 1 J_ A A/1 

— U.l i it U.U4 


A 1 1 J- A A/1 

U.il ± U.U4 




■3 


1 J.J ± Z.J 


A AO J- A A/I 

—U.UZ it U.U4 


A 1 C 1 A A1 

U.l J it U.Ui 


A 1 ■I J_ A A/1 

—U.l J it U.U4 


A AO J- A A/1 

U.UZ it U.U4 




4 


1 a 1 4- 7 
1 J.l ±A.I 


A 1 A _|_ A A/1 
U. lU it U.U4 


n 00 -1- n no 

U.ZZ It U.UZ 


_n in -U n c\a 

— U.IU it U.U4 


n 1 -1- n (\A 

U. iZ it U.U4 


NGC 4649 


1 


18.2 ±2.5 


0.04 ± 0.04 


0.28 ± 0.02 


-0.22 ±0.04 


0.06 ± 0.04 




z 


10 4-0 S 

ly.z It z.o 


U.UU It u.u^ 


n 00 4- n no 
u.zy It U.UZ 


_n 00 4- n (\A 

—U.ZZ It U.U4 


n no 4- n nA 
U.U / It u.u^ 




■J 


of\ n -1- < 
ZU.U it Z.J 


A no 4_ A no 

— U.Uo it U.Uj 


f\ OA 1 n m 
U.ZU it U.Ui 


n 0'2 4- O C\1 
— U.Zj It U.Uj 


n AO 4_ n AO 
—U.UJ it U.UJ 




4 


1 Q "2 -1- Q 
iO.J it Z.O 


U.Uj it U.U4 


O OO 4- O f\0 

U.zy it U.UZ 


n OO 4- O f\A 
—U.ZZ it U.U4 


A OO 4- A 0/1 
U.U / It U.U4 


NGC 4697 


1 


15.6 ±2.2 


-0.29 ±0.03 


0.08 ±0.03 


-0.36 ±0.04 


-0.28 ±0.03 




2 


15.5 ±2.2 


-0.29 ±0.03 


0.08 ±0.03 


-0.35 ±0.04 


-0.27 ±0.03 




3 


15.6±2.1 


-0.31 ±0.04 


0.06 ±0.02 


-0.36 ±0.04 


-0.30 ±0.04 




4 


15.6±2.5 


-0.29 ±0.04 


0.09 ±0.03 


-0.37 ±0.05 


-0.28 ±0.04 


NGC 5638 


1 


11.6±2.0 


-0.05 ±0.03 


0.13 ±0.02 


-0.17 ±0.04 


-0.04 ±0.03 




2 


11.7±2.1 


-0.06 ±0.03 


0.13 ±0.02 


-0.16 ±0.03 


-0.03 ±0.03 




3 


12.1 ±2.5 


-0.10 ±0.03 


0.09 ±0.01 


-0.17 ±0.03 


-0.08 ±0.03 




4 


11.6±2.1 


-0.04 ±0.04 


0.13 ±0.02 


-0.16 ±0.04 


-0.03 ±0.04 


NGC 5812 


1 


7.2±1.4 


0.22 ±0.04 


0.16 ±0.02 


0.07 ±0.04 


0.23 ±0.04 




2 


7.5±1.1 


0.20 ±0.03 


0.16 ±0.02 


0.08 ± 0.03 


0.24 ±0.03 




3 


8.2±1.3 


0.14 ±0.04 


0.11 ±0.01 


0.06 ± 0.04 


0.17 ±0.04 




4 


7.0±1.3 


0.23 ±0.04 


0.17 ±0.02 


0.07 ±0.04 


0.24 ± 0.04 


NGC 5813 


1 


24.3 ±2.1 


-0.21 ±0.05 


0.12 ±0.03 


-0.32 ± 0.06 


-0.20 ± 0.05 




2 


24.4 ±2.0 


-0.22 ±0.05 


0.12 ±0.03 


-0.31 ±0.06 


-0.19 ±0.05 




3 


24.7 ±2.0 


-0.27 ±0.05 


0.09 ±0.02 


-0.34 ±0.05 


-0.25 ±0.05 




4 


24.3 ±2.0 


-0.20 ±0.06 


0.12 ±0.03 


-0.31 ±0.07 


-0.19 ±0.06 



TRACER ETAL. 



Table 6b — Continued 



Name 


Model 


Age (Gyr) 


[Z/H] 


[E/Fe] 


[Fe/H] 


[E/H] 




1 

1 




n 1 7 -u n 07 


n 1 9 -1- n 09 

L/. iZ m L/.UZ 


Oft -t- 07 
U.L/D HI U.L// 


1 R -1- 07 
L/. io in L/.U / 




2 


4 -t- 1 9 


1 ^ 4- 07 


n 1 9 4- n 09 

w. iz m w.uz 


Ofi -t- 07 

u.uD in u.u/ 


1 R -}- 07 




3 


4.8 ± 1.3 


0.11 ±0.04 


0.08 ±0.01 


0.05 ±0.04 


0.13 ±0.04 




4 


4.2±1.0 


0.18 ±0.06 


0.13 ±0.02 


0.06 ±0.06 


0.19 ±0.06 




[ 


23 1 i 2 5 


— V^. 1^ ZLL V^.UH- 


n 1 Q -1- n 03 


_0 1 -t- OS 
— u.ji in u.uj 


_0 1 9 -1- 04 
u. iz m u.uH 




Z 


9'^ A -1- 9 zl 

Z.).^ Ill Z.4 


1 -t- 04 
U. ID III U.U't 


1 Q -t- 0'^ 


1 -t- OS 


1 9 -t- 04 
U. IZ III U.U^ 




3 


23.6 ±2.3 


-0.21 ±0.06 


0.13 ±0.02 
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Correcled values: 
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Fig. 3. — Line strengths of ear ly-typ e galaxies in the Gonzalez (1993) sample in the central Tc/S aperture, corrected to solar abundance ratios ([E/Fe] = OXusing 
the method described in Section 3.1.3 and enhancement model 4. Solar abundance model grids from Worthey (1994) are again superimposed as in Figure 1. The 
corrections to solar abundance ratios bring inferred SSP ages and metallicities into good agreement between the Mgfc-H/3 and the (Fe)-H/3 diagrams, as expected. 
These figures indicate the final central SSP values for t and [Z/H]. 



Corrected values: r^/2 




0,4 0.5 O.B 0.7 0,9 0.^ 0,5 0,6 



log Mgb log <Fc> 



Fig. 4. — Line strengths of early-tynp galaxies in the Gonzalez (1993) sample in the global re {2 aperture, corrected to solar abundance jatios ([E/Fe] = 0) using 
the method described in Section 3.1.2 and enhancement model 4. Model grids from Worthey (1994) are again superimposed as in Figure hi Ages and metallicities 
in the two panels agree, as they do tor the central values in Figure H (with the exception of NGC 720, which lies far off the grid). These ligures indicate the final 
global SSP values for t and [Z/H]. 
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Fig. 5. — SSP ages, metallicities, and enhancement ratios in models 1, 2, and 3 compared to model 4 (which has C+0+). (tlata are for the central r^/S aperture.) 
Models 1, 2, and 4 give nearly identical population parameters because the C abundance changes little among them (cf. Table W); model 3 deviates strongly from the 
others (especially in age and enhancement) because of the strong dependence of MgZ) on C abundance (TB95; cf. Tablen). 
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Fig. 6. — A comparison of measured and model-derived indices for the four element enhancement models. Shown here are three new C-sensitive Une indices for 
the central regions (from Lick/IDS data), plus Mgb. shown for reference. The three new indices were not used to compute the original SSP parameters and therefore 
provide an independent check on them; "computed" indices are those predicted by the SSP parameters. Model 3 (open squares) systematically underestimates the 
C-sensitive indices compai'ed to models 1, 2, and 4. This is because its C abundance is too low. The good agreement for Mgb shows that we have the proper 
solutions. 
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4. SSP-EQUIVALENT PARAMETERS FOR THE G93 SAMPLE 

This section presents a brief overview of the resultant SSP- 
equivalent population parameters for the G93 galaxies; detailed 
discussion is reserved to Papers II and III. Our focus here is 
on the preferred model 4 (C and O both up), but results from 
models 1 and 2 are similar (model 3 being ruled out). 

Figure presents histograms of f, [Z/H], and [E/Fe] for the 
G93 sample through the r^/S aperture. The original conclu- 
sions of G93 are confirmed using this more rigorous analysis: 
the central stellar populations of galaxies in this sample span a 
large range of SSP-equivalent ages, from 1.5 f (Gyr) 18 
(more than 1 dex), but a relatively narrow range in [Z/H], 
-0.1 < [Z/H] ^ +0.6, and an even smaller spread in [E/Fe]. 
The metallicity distribution has a peak at ([Z/H]} = +0.24 and 
a dispersion of <t{[Z/H]) = 0.14, while the enhancement dis- 
tribution peaks strongly at ([E/Fe]) = +0.20 with a dispersion 
cr([E/Fe]) of only 0.05 (these values vary slightly with the 
model). 

A striking fact to emerge from Figure^ is how mild the mean 
metallicities and enhancements of ellipticals really are. Match- 
ing the high Mg i ndex values of elHpticals has been problematic 
in the past (e.g., Matteucci 1994; Greggio 1997), and previ- 
ous authors have typically invoked rather large enhancements 
in th e range [E/Fe] = +0.3-0.5 (W eiss, Peletier & Matteucci 
1995; [Trager 1997| ; Preggio 1997| ). With the TB95 response 
functions, however, the average [Z/H] is only a factor of two 
higher than solar, and the average [E/Fe] is only +0.2. The lat- 
ter is small compared to the maximum value of [a/Fe] ^ +0.5 
foun d in metal-poor Galactic stars (W heeler, Sneden & Truran 
1989; pSdvardsson et al. 1993| ), which is widely regarded as an 
empirical upper limit to the amount of depression in Fe that 
can result from total suppression of SNae la. The depression 
of the Fe-peak in ellipticals appears to be much less than this 
and should be easier to accommodate with reasonable galacto- 
nucleosynthesis models. 

Figure H presents similar histograms of f, [Z/H], and [E/Fe] 
for the aperture. The global stellar populations span a 
slightly larger range of ages, from 1.5 to 25 Gyr and a slightly 
wider range of metallicities, -0.3 <; [Z/H] ^ +0.7 (with NGC 
720 at [Z/H] = 1.05), although some of this larger scatter is 
surely due to the larger uncertainties in the rp/2 line strengths. 
Otherwise, the shapes of the distributions are similar Com- 
paring r(,/8 with rp/2 shows that mean ([Z/H]) is down by 
0.18 dex in the outer parts, indicating that the outer regions are 
slightly more metal-poor than the centers. The outer mean en- 
hancement ([E/Fe]) is lower by only 0.03 dex, however, con- 
firming the conclusion of Worthey et al. (1992), Davies, Sadler 
& Peletier (1993), and G93 that enhancement gradients within 
galaxies are weak. Ages increase slightly outwards, the outer 
parts being on average roughly 25% older Overall, the differ- 
ences among galaxies are much more striking than the differ- 
ences within galaxies, at least in the G93 sample, through these 
apertures. 

5. UNCERTAINTffiS AND SYSTEMATIC ERRORS 

This section assesses both zeropoint and scale errors in f, 
[Z/H], and [E/Fe]. We begin by examining our basic as- 
sumption that the ages, metallicities, and enhancement ratios 
we have derived above represent true light-weighted ages and 
abundances of elliptical galaxies. In particular, we first ask 



whether the apparent large age spread among the G93 galax- 
ies could be due to spurious effects. 

5.1. Hf] as an age indicator 

The assumption that we are measuring real ages of stellar 
populations rests on the further assumption that H/3 light is 
coming purely from main-sequence and red giant-branch stars. 
We now discuss three scenarios whereby H/3 might be contam- 
inated by light from other sources. 

(1) Fill-in by emission (see Section 2.2.2). The extreme 
form of this hypothesis says that all ellipticals are actually 
young and that the apparent large age spread is due entirely 
to variable amounts of infill by emission. This extreme view 
is strictly ruled out by numerous observational studies of emis- 
sion in elliptical galaxies. For example, G93's plot of preci- 
sion continuum-subtracted spectra (G93, Figure 4.10) shows 
that emission is nearly always less than a few tenths of an A, 
not nearly large enough to create the observed age spread. In 
the same vein, Carrasco et al. (1996) went so far as to suggest 
that no emission corrections should be applied at all to most 
ellipticals, implying that any emission can at most be small. 
The final poi nt is that H/3 co rrelates strongly both with Mg2 
and a (G93; J0rgensen 1997), inconsistent with emission fill- 
in, which varies unpredictably from galaxy to galaxy. 

A more reasonable hypothesis is that errors in the emission 
correction contribute noticea bly to the age spread. Such errors 
were investigated in Section 2.2.2, where we noted that scatter 
in the H/3/[0 III] ratio would induce age errors of only ±9% 
for typicalgalaxies. An even more drastic test is presented in 
Appendix pi which shows that neglecting the emission correc- 
tion altogether affects a few strong- [O III] galaxies but makes 
at most small changes in the broad age distribution. 

(2) Contamination by blue horizontal branch stars (BHBs). 
BHB stars are not present in the standard Worthey (1994) mod- 
els, which assume red clumps for old metal-rich populations. 
BHB stars might come from an anomalous BHB population as- 
sociated with the metal-rich stars, or from contamination by a 
normal BHB associated with a subordinate metal-poor popu- 
lation. By "BHB," we mean blue horizontal branches similar 
to M 92, which would contribute significantly to the light at 
4000-5500 A, not the extremely hot horizontal branches iden- 
tified in populations like NGC 6791 (Liebert, Saffer & Green 
1994) which contribute primarily to 1 500 A flux and the "UV 



upturn" (e.g., Lee 1994 



Yietal. 1999). 



The galaxy M 32 can be used to rule out the hypothesis that 
BHBs alone are responsible for the large H/3 excesses seen 
in high-Hp ellipticals. It can be shown that nearly the entire 
red clump in M 32 would have to be moved to a BHB at ap- 
proximately spectral type mid-F to explain its high H/3 index 
( Burstein et al. 1984 ); this is strictly ruled out by blue spec- 
ti'al indices ( |Rose 1985 , 1994). Moreover, the HB has actually 
been de tected in the outer part of M 32 by HST (G rillmair et 
al. 1996) and is seen to be mostly red.' Extrapolating the G93 
indices outward to this field and matching to W94 models there 
yields an excellent fit to b oth the integrated col ors and the color 
of the RGB at this point (Grillmair et al. 1990, supporting the 
assumption that the HB is indeed red. 

The existence of a dominant BHB population in metal-rich 
ellipticals is not expected on astrophysical grounds. If ellipti- 
cals were very old (i.e., > 18 Gyr, as Lee 1994 has suggested), 

'To be precise, the Grillmair et al. (1996) data are not deep enough to rule out a small number of BHB stars ( [jrillmair et al. 199^ ; C. Gallart, priv. comm.). but 
the lack of point sources in archival F300W images suggests that any BHB must indeed be weak. 
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Fig. 7. — Distributions of central (r^/S) stellar population parameters for the Gonzalez (1993) sample using enhancement model 4. As expected from Figures 1 
and 3, these local field ellipticals span a wide range of ages but a small range (< 0.3 dex) of metallicities and a very small range in abundance enhancements. 
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Fig. 8. — Distributions of global (re/2) stellar population parameters for the Gonzalez (1993) sample using enhancement model 4. The global parameters here are 
slightly offset to older ages and lower metallicities from the central parameters in Figure 7, but otherwise the distributions are strikingly similar. 



then BHBs could conceivably be significant components, but 
such large ages viol ate current constrain ts on the age of the 
Universe (see, e.g., aratton et al. 1997). No solar metallic- 



ity cl uster populations in the Milky Way have BHBs (W orthey 
1994), although we note that Rich et al. (1997) have discov- 
ered significant M3-like BHB populations in two metal-rich 
Galactic bulge globular clusters, NGC 6388 and NGC 6441 
([Fe/H] ^ -0.5). However, these two globulars are the dens- 
est known in the Galactic globular cluster system; th e fact that 
BHB stars occur precisely there (Sosin et al. 1997) suggests 
that dynamical interactions are the cause. We conclude that the 
occurrence of BHB stars in low-density systems like giant el- 
liptical galaxies is unlikely, but a deeper understanding of their 
presence in these globulars is obviously necessary. 

Contamination by BHBs from a subordinate metal-poor pop- 
ulation also does not seem probable. As noted, such contami- 
nation by a trace blue BHB component cannot materially affect 
the indices of high-H(3 galaxies, but perturbations in weak-Hf] 
galaxies should be considered. For example, 5% of the V-band 
light in metal-poor BHB stars would decrease the inferred age 
of a galaxy from 13 Gyr to 8 Gyr at solar metallicity. However, 
Rose (1985, 1994), using a set of high-resolution spectral in- 
dices in the 4000 A region, has shown in M 32 and eight strong- 
lined ellipticals that no more than ~ 5% of the light in the blue 
region (and less than 2% in the V-band) can come from very hot 
stars (FO and earlier). This falls short by a factor of two. More- 
over 5% of V-band light in BHB stars would imply that alto- 
gether ^25% of the total light would have to come from metal- 
poor stars. This is twenty-five times more than the amount of 
metal-poor ([Z/H] < -1.5) V-band light actually found in the 
outer part of M 32 by Grillmair et al. (1996). That a much larger 
quantity of metal-poor stars could be found near the centers of 
more metal-rich elliptical galaxies seems implausible. 

We conclude that contamination by metal-poor populations 



is a negligible perturbation to the central ages of the G93 galax- 
ies and could cause at most a ^10% reduction near r^, if that. 

(3) Contamination by blue stragg ler stars (BS Ss). A typical 
BSS has My « 3 mag, B- V « 0.2 ( $ailyn 1995] ), and specti-al 
type A8-F0. From Worthey et al. (1994), dwarf A8-F0 stars 
have up « 5.8-7 A. To explain the high H/3 sti-ength of M 32 
(H/3 = 2.4 A) as arising from a population of blue stragglers 
superimposed on an old (15 Gyr), solar- metallicity population 
(H/3 = 1.5 A) would require that ^ 15% of the V-band Hght 
come from BSSs. This implies a BSS specific frequency of 
^ 275 per 10** Lq, which is a factor of 8 higher than seen in the 
most BSS-rich Galactic globular cluster (Palomar 5) and a fac- 
tor of 28 higher than seen in the average G alactic globular clus- 
ter ( l^erraro, Fusi Pecci & Bellazzini 1995 ). We again conclude 
that high-Hj3 galaxies like M 32 are immune to perturbations 
by spurious hot components such as BSS stars. 

Consider next a trace contamination by BSSs in low-H/3 
galaxies. For example, to decrease the age of an elliptical from 
13 Gyr to 8 Gyr at solar metallicity, BSSs would again need 
to contribute ^ 5% of the V-band light. For NGC 3379, this 
implies a specific frequency of ^ 120 per lO** L0. This spe- 
cific frequency is a little more than a factor of 3 higher than that 
seen in the most BSS-rich G alactic globular cluster (F erraro, 
Fusi Pecci & Bellazzini 1995). In the absence of a complete 
theory of BSS formation, a factor of three increase might not 
be impossible. On the other hand, as noted. Rose (1985, 1994) 
has shown that no more than ~ 2% of V-band light can come 
from hot stars FO and earlier in M 32 and eight strong-lined 
ellipticals. This is less than half the light required and would 
perturb the age from 1 3 Gyr to only 1 1 Gyr, a reduction of only 
15%. 

To summarize, the leading hot-star contaminants, BHBs and 
BSSs, both peak at temperatures at or hotter than FO, whereas 
the blue line-strength data of Rose (1985, 1994) imply that the 
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great bulk of Balmer absorption must be coming from cooler 
F and G stars, at or at most only slightly hotter than the de- 
rived turnoff temperatures. Barring some as-yet-undiscovered 
contaminating population of cooler stars, the Rose limits imply 
that contamination of H/3 by non-main sequence stars can re- 
duce the ages of even the oldest ellipticals by at most 10-15%. 

5.2. Errors due to theoretical model uncertainties 

The next three sections assess additional sources of system- 
atic errors; results are collected in Table ^ This section dis- 
cusses theoretical model uncertainties caused by errors in the 
stellar isochrones and line-strength response functions of TB95. 
The major uncertainty in the interior models is the age scale, 
which is continually being refined. With the recent release of 
parallaxes from the HIPPARCOS satellite, much effort has been 
spent recalibrating the ages of Galactic globular clusters using 
both the ne w parallaxe s and up-to-date models o f stellar evo- 
lution ( e.g.,|Reid 199% 1998; Pratton et al. 1997|; C haboyer et 
al. 1 998; IGrundahl, Vandenberg & Andersen 1998t P ont et al. 
1998; ISalaris & Weiss 199S| ). This effort has brought the ages 
of the oldest globular clusters down from ~ 14-15 Gyr to ~ 12 
Gyr, a reduction of 15%. At least half of this reduction is due 
to corrections in the metallicity scale of globular clusters and to 
the use of more up-to-date stellar evolutionary models (G ratton 
et al. 1997). 

These age redeterminations have so far been restricted to 
clusters with metallicities [Fe/H] <^ -0.7. At the metallicities 
typical of elliptical galaxies, the effect of the age recalibra- 
tions is not yet known but could be as much as ^ 20%, just 
by using isochrones from the most modern stellar evolutionary 
models. This agrees with Chariot, Worthey & Bressan (1996), 
who found that absolute ages are uncertain at the 25% level 
in stellar populations with ages > 10 Gyr, resulting almost en- 
tirely from the choice of different stellar models. Below and 
in Appendix ^ we explore the effect of substituting "Padova" 
isochrones by Bertelli et al. (1994) for those of W94 and find 
that young ages differ by 35% but that old ages change by only 
4%. As a rough rule of thumb, we assume that both the age 
zeropoint and age scale of the models are uncertain at the ^ 20- 
25% level. 

The effect of errors in the theoretical response functions of 
TB95 is illustrated in Figure ^ which is a schematic repeat of 
Figure 1(b) showing H/3 versus (Fe). Figure^ shows a galaxy 
plotted two ways, one using raw (Fe), the other using the value 
of (Fe) inferred from Mgb by assuming solar abundance ratios 
(call this (Fe(Mg))). (Fe) Hes to the left of (Fe(Mg)), indi- 
cating Fe depression. Applying the TB95 corrections for non- 
solar [E/Fe] moves (Fe) to the right and (Fe(Mg))to the left, 
as shown by the arrows (the correction to H/3 is small and is 
ignored). When the correct value of [E/Fe] is reached, the two 
points coincide, giving final f , Z, and [E/Fe] (right hand panel). 

Where the solution lands is evidently governed by the rel- 
ative lengths of the two correction vectors; for model 4, this 
ratio is Alog(Fe(Mg))/Alog(Fe) = 1.25. The systematic er- 
rors of the final point depend mainly on the error of this ratio. 
Assuming that the two response functions of TB95 are individ- 
ually uncertain by as much as 30% and that the errors of their 
three stellar types add in quadrature, the resultant zeropoint un- 
certainties are 3% in age, 0.10 dex in [Z/H], and 0.04 dex in 
[E/Fe] for highly enhanced galaxies. Since this last error drops 
to zero for galaxies with [E/Fe] = 0, we derive an overall scale 
uncertainty in [E/Fe] of <20%. 

We note that fundamental uncertainties in stellar models, for 



example the use of a single-parameter mixing length theory for 
convection or the detailed effects of rotation and diffusion, may 
induce additional, unknown systematic errors in our absolute 
age estimates. Such uncertainties also affect the globular clus- 
ter age scale. At present, our estimated uncertainties in the ab- 
solute ages of galaxies therefore should be considered to be rel- 
ative to the globular cluster age scale. 

5.3. Errors due to empirical model uncertainties 

Errors in this category include errors in the metallicity and 
temperature scales of the Lick/IDS fitting functions and errors 
in the fitting accuracy of the functions themselves. 

We have checked the metallicity scale of the Lick/IDS sys- 
tem by comparing our assum ed stellar [Z/H] values (summa- 
rized by Worthey et al. 1994 ) with the compilation of published 
spectroscopically-determined values by Cayrel de Strobel et al. 
(1997). For stars with log^ > 4 (mostly dwarfs), the Lick/IDS 
metallicity scale is in excellent agreement with the published 
values at all [Fe/H]. For giants with [Fe/H] <^ 0, the Lick/IDS 
metallicity scale is within 0.05-0.1 dex of the Cayrel de Strobel 
scale (systematically slightly high). However, for giants with 
[Fe/H] > ("SMR" stars), the Lick/IDS metallicity scale de- 
viates strongly from the Cayrel de Strobel scale, such that the 
Lick/IDS giants appear to be more metal rich. The Lick/IDS 
metallicity scale for giants is based on the narrow-band photo- 
metric metallicity scales of Hansen & Kjsergaard (1971) and 
Gottlieb & Bell (1971), and on the high-resolution spectral 
study o f Gustafsson, Kjsergaard & Anderson (1974) ( Faber et 
al. 1985). In contrast, the Cayrel de Strobel et al. (1997) cata- 
log is populated in the SMR giant regime by older spectroscopic 
abundance determinations based on stellar atmospheres that 
typic ally (1) have a too-low solar iron abundance (M c William 
1997) and (2) d o not properly acc ount for molecule formation 
in SMR giants (Castro et al. 1996). Correcting the abundances 
of SMR giants for these two effects suggests that the Lick/IDS 
scale may actually be very close (possibly 0.05-0.1 dex too 
high) to the modern spectrosco pic metallicity scale, even at 
[Fe/H] 0.4 ( |Castro et al. 1996| ; Mc William, priv. comm.). 

The next question is whether the fitting functions are in fact 
good fits to the stellar line strengths. By inspecting the resid- 
ual diagrams in Gorgas et al. (1993) and Worthey et al. (1994), 
we estimate that any systematic errors in the metal-line fits are 
less than 3%, which translates to zeropoint uncertainties of 0.05 
in [Z/H] and 0.10 in [E/Fe] (the former averages Fe and Mg 
while the latter differences them, accounting for its larger er- 
ror). The crucial function for age is the fit to H/3 versus V -K 
for main sequence A-F stars. Again, we estimate that the ba- 
sic line-strength calibration level is accurate to better than 3% 
in this interval, which translates to about 10% in age. Finally, 
the temperature scale {V -K vs. Te) of main sequence stars is 
needed to attach H/3 strengths to the t heoretical isochrones. An 
error of 100 K (Worthey et al. 1994) again translates to about 
10% in age. Note that all these errors in the fitting functions af- 
fect only the absolute zeropoints of age, metallicity, and [E/Fe] 
but not their differential values. 

5.4. Errors due to unknown element enhancements 

Our treatment of element enhancements is crude — we simply 
group all elements into three categories (enhanced, depressed, 
and fixed) and assume that differences within each group are 
nil. The group assignment of certain elements is also uncertain. 
Unknown element abundance ratios introduce errors in the pre- 
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Fig. 9. — A schematic representation showing how SSP pai'ameters are derived for a galaxy with non-solar abundance ratios. A model non-solar galaxy with 
(t, [Z/H], [E/Fe]) = (12,-0.25,0.3) is plotted in panel (a) with its true value of (Fe) (solid square) and the value of (Fe) that would be inferred from its Mgb strength 
assuming solar abundance ratios (solid triangle; called (Fe(Mg)) in the text). This latter point lies to the right of (Fe) reflecting non-solar ratios. The vectors show 
corrections to each feature for the non-solar enhancement ([E/Fe] ^0; small changes in H/3 are ignored). When the correct enhancement is chosen, the two vector 
tips agree. The sum of the two vectors determines the size of the enhancement, while their ratio determines the location of the final point, which gives t and [Z/H]. 
The corrected point is also shown in panel (b) and coincides with the inferred (f,[Z/H]) of (12,-0.25). Errors in the two coiTection vectors affect the final population 
parameters as described in the text. 



dieted index response functions and, to a smaller extent, in the 
theoretical stellar evolutionary tracks. 

According to TB95, the elements that significantly influence 
the indices used here are Fe, Cr, C, and Mg. Fe and Cr are pro- 
duced in both Type la and in interm ediate-mass progenitor Type 
II SNae (Woosley & Weaver 1995). They should vary closely 



together by virtue of similar nucleosynthesis; i.e., their relative 
uncertainty should be small. Breaking the link between Fe and 
Cr, for example by decreasing [Cr/Fe], would have the effect of 
altering the Mgb index strength without significantly affecting 
other indices (TB95). In our own galaxy, however, [Cr/Fe] i s 
solidly at the solar value until [Fe/H] -2 ( fvlcWilHam 1997[ ), 
much lower than the metallicities of interest in elliptical galax- 
ies. We will discuss possible element-to-element variations in a 
future paper. Likewise we have tested the sensitivity of the in- 
dices to C explicitly in models 1-4 and found that low-C models 
(like model 3) are ruled out. With this eliminated, remaining 
uncertainties due to the C abundance variations are limited to 
10% in age, 0.05 dex in [Z/H], and 0.01 dex in [E/Fe] (Fig- 
ure |). 

A larger source of uncertainty arises from uncertain ratios 
within the Type II SNae group. The metallicity [Z/H] is con- 
trolled by O, which has little spectroscopic signature (TB95), 
while a major spectral impact comes from Mg. Our inferred 
values of [Z/H] thus depend critically on the assumption that 
Mg and O track one another. Breaking this link, e.g., by en- 
hancing Mg over O, could reduce our inferred [Z/H]'s sub- 
stantially. For example, suppose that [O/H] is always solar re- 
gardless of Mg (this would place O in the depressed group in 
metal-rich galaxies). Since O contributes half the mass in Z 
(see Table Bh, our values of [Z/H] would be overestimated by 
a factor of two. Correcting for this would reduce the typical 
[Z/H] from 0.26 to 0.13, and in so doing would increase ages 



by about 20%; enhancements [E/Fe] would remain unchanged. 
We are thus relying quite heavily on the notion that decoupling 
O and Mg is astrophysically unreasonable. 

Finally, our analysis assumes that isochrone shape and loca- 
tion are unaffected by the exact value of [E/Fe] or by the de- 
tailed pattern of element enhancements within [E/Fe]. Existing 
models suggest that this assumption might be tolerable. Salaris 
& Weiss (1998) have calculated an isochrone for an old popu- 
lation model with [Z/H] = -0.3, [E/Fe] = +0.4, and non-solar 
['^hpe/'^lpe] = 0.12. Log at the turnoff shifts to the blue by 
0.0044, while log 7; on the RGB shifts to the blue by 0.01 1 rel- 
ative to a scaled solar model. The blueward shifts should scale 
in proportion to both [Z/H] and [E/Fe], while the shape change 
may also depend on [Xhpe/^lpe] (S alaris, Chieffi & Straniero 
1993; [Salaris & Weiss 1998| ). A typical G93 galaxy is four 
times more metal-rich than their model but smaller by a factor 
of two in [E/Fe]. The quantity [^hpe/^lpe] is also likely to 
be smaller, being -1-0.07 in model 4 (Table ^ versus -1-0.12 in 
their model. On balance, the net shifts and shape changes in the 
elliptical isochrones are plausibly no more than twice those in 
their model. 

The effects of such motions would be small. A shift of log 
Tf = 0.0044 at the turnoff causes a change of only 0.016 in log 
H/3, for a change in age of 6%. A shift of log Tg = 0.01 1 on the 
RGB causes a change in metal lines of the same amount, for a 
change in [Z/H] of about 0.05 and no change in [E/Fe]. Even if 
multiplied by two, as estimated above, these effects would still 
be small compared to other errors. On tho other hand, it should 
be stressed that the effect of non-solar ratios on isochrone lo- 
cation is acclerating at high metallicity, and the above models 
were calculated for metallicities considerably smaller than what 
we require. A failure of O to track Mg (as mentioned above) 
could also introduce further shape changes that have not yet 
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been modeled in detail. In sum, our assumption that isochrone 
shape and location are unaffected by the value of [E/Fe] or 
by the pattern of non-solar enhancements within [E/Fe] looks 
promising but is in need of further validation. 

5.5. Error summary 

The results of the preceding sections, plus some additional 
experiments in Appendices ^ and ^, are summarized in Ta- 
ble Age errors are significant — several terms amount indi- 
vidually to 10-25% and their addition is uncertain. Some age 
errors are also larger for weak-H/3 objects and therefore tend to 
stretch or compress the age scale. However, most of the errors, 
including those in age, are simple zeropoint shifts. Future ap- 
plications will take advantage of the relative robustness and use 
the data differentially. 

The galaxy M 32 offers a final check on the zeropoints of 
both [Z/H] and [E/Fe]. The integrated spectrum of M 32 
has been modeled by many authors, and the upper CM dia- 
gra m of the outer parts has been measured (G rillmair et al. 
1996). All spectrum modelers concur that a mix of moderately 
young stars of near-solar metallicity matches every known fea- 
ture of the spectrum. The mean turnoff spectral type within 



the re/2 aperture is accuratel y known to be F7-8 (Faber 1972 
O'Connell 1980| ; ^ose 1994| ), while the light-weighted metal- 
licity in the Grillmair field is [Fe/H] = -0.25. The enhance- 
ment ratio [E/Fe] is also known to be small based on the ex- 
celle nt spectral fi ts using solar-neighborhood-abundance stars 
(e.g., [Faber 197^). 



These independently measured parameters agree well with 
the SSP-equivalent parameters. G93 indices for the re/2 aper- 
ture yield an SSP-equivalent age of 5 Gyr, a mean [Z/H] of 
-0.07, and an [E/Fe] of -0.07. These parameters imply a 
turnoff spectral type of exactly F7-8, as the modelers have con- 
cluded, and the near-solar [Z/H] and [E/Fe] also agree with 
their results. Extrapolating the G93 indices outward, Grillmair 
et al. (1996) found an SSP-equivalent age in their field of 8 Gyr, 
a mean [Z/H] of -0.25, and [E/Fe] of -0.05. This metallicity 
coincides precisely with the metallicity distribution that they 
inferred from the color locus of the RGB for that assumed age. 
Putting this information together, we conclude that the actual 
absolute uncertainties in both [Z/H] and [E/Fe] are < 0.05, at 
least for galaxies close to solar composition like M 32. 

6. COMPARISONS WITH PREVIOUS STUDIES: EVIDENCE FOR 
INTERMEDIATE- AGE POPULATIONS 

Many previous studies have examined the line strengths and 
colors of elliptical galaxies to determine their stellar content. 
A complete review of all previous models is beyond the scope 
of this paper (see Chariot, Worthey & Bressan 1996, Vazdekis 
et al. 1996, and Arimoto 1996 for comparison of some mod- 
ern stellar population synthesis models, and Worthey 1998 for 
a historical review of the metallicities and abundance ratios of 
early-type galaxies). We concentrate here on previous investi- 
gations that derived SSP ages and models by using the Balmer 
lines. We begin with the results of TCB98 and then turn to those 
of other workers. Consideration of other methods, in particular 
those using colors, is delayed to future papers. 

6.1. Model dependence of derived stellar population 
parameters: comparison with TCB98 



In a recent paper, TCB98 have analyzed line strengths of 
the G93 galaxies in the context of their own stellar population 
models. These models are based on isochrones by Bertelli et 
al. (1994) (which, like ours, neglect the effects of [E/Fe] j^O), 
the original fitting functions for H/3 and (Fe) by Worthey et al. 
(1994) (which also neglect [E/Fe] ^0), and a fitting function 
for Mg2 by Borges et al. (1995), who claim to take into account 
the effects of [E/Fe] y^Q. Like us, TCB98 assume that their 
isochrones depend only on bulk metallicity [Z/H] but not on 
[E/Fe] (which they call [a/Fe]). To determine the stellar pop- 
ulation parameters logf, [Z/H], and [a/Fe], TCB98 compute 
averaged derivatives (from their models) of (the logarithms of) 
Mg2, (Fe), and H/3 versus population parameters. These deriva- 
tives are then inverted to derive a series of linear equations that 
yield relative values of logf, [Z/H], and [a/Fe] as functions of 
Mg2, (Fe), and H/3. This solution method is equivalent to as- 
suming that all line strengths depend linearly on all parameters, 
which is marginally inconsistent with the curved shapes of the 
actual grids (cf. Figure [l|). 

Figure [l^ (top row) shows the stellar population parameters 
derived by TCB98 as a function of our derived parameters. For 
this comparison we have used our enhancement model 3 with 
C depressed and O enhanced, which is closest to their model. 
While the two studies roughly agree in the distribution of ages 
of the G93 galaxies, they are discordant in both [Z/H] and 
[E/Fe], for which the slopes of their values versus ours deviate 
strongly from unity. These different inferred [Z/H] and [E/Fe] 
values imply different interpretations of the star formation his- 
tories of these galaxies (particularly in the logf-[Z/H] relation; 
see below). 

To isolate the source of the differences, we constructed new 
models ("Padova") by substituting the Bertelli et al. (1994) 
isochrones for the RYI/VandenBerg isochrones used by W94 



(Appendix JA). The results are presented in the bottom row 
of Figure |^ The match between the new Padova models 
and the standard W94 models is quite good, with only slight 
slope changes and mild zeropoint offsets; logf decreases by 
from 10% to 30%, [Z/H] increases by less than +0.08 dex, and 
[E/Fe] increases by no more than +0.02 dex in the Padova mod- 
els. 

The differences between the present results and those of 
TCB98 in the top row are therefore not caused by differences 
in the isochrones, but must rather stem from one of the follow- 
ing other differences: (1) use of different response functions 
for [E/Fe] ^ for Mg2, Fe5270, and Fe5335; (2) use of Mg2 in- 
stead of Mgb; and/or (3) use of a linearized solution method for 
deriving ages, metallicities, and enhancement ratios from ob- 
served line strengths. Inspection shows that differences (2) and 
(3) are most likely minor; in particular, the G93 Mg2 strengths 
are less reliable than the Mgb strengths, but broadly the two 
agree fairly well. Likewise, the linearized method deviates at 
large distances from the middles of the grids owing to grid cur- 
vature, but these differences are not large enough to cause the 
global slope differences seen in Figure HQ. 

Difference (1), the use of different response functions, dom- 
inates the differences in [Z/H] and [E/Fe]. TCB98 use the 
Borges et al. (1995) fitting function for Mg2, which nominally 
takes [E/Fe] ^ into account but is derived from only a small 
set of calibration stars. Furthermore, by using the original 
Worthey et al. (1994) fitting functions for Fe5270 and Fe5335 
without correction for [E/Fe] j^Q, they implicitly assume that 

^The stellar metallicil£s_are most likely on a different metallicity scale, as they are drawn directly from Cayrel de Strobel et al. (1997), which may be unreliable at 
[Fe/H] > 0; see Section p.4[ The Borges et al. sample is also deficient in calibrating stars on the RGB. 
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Table 7 

Summary of Systematic Errors 



Source 


Amount 


Effect 








Metal-poor BHB contamination 


< -10%^ 


Weak-H/3 objects only 


Blue straggler contamination 


<-15% 


Weak-H/3 objects only 


Theoretical isochrones 


20-25% 


Zeropoint 


Theoretical isochrones 


20-25% 


Scale 


Theoretical response functions (TB95) 


3% 


Zeropoint 


Up fitting function vs. V-^T 


10% 


Mainly zeropoint 


V — K vs. Te 


10% 




Unknown O/Mg ratio 


20%'^ 


Metal-rich objects only 


Unknown C abundance 


10% 


Zeropoint 


Effect of [E/Fe] 7;'0 on isochrones 


6-12% 


Zeropoint 


Overestimation of H/3 emission correction 


H-3% 


Zeropoint 


Undercorrected emission fill-in 


< +25% 


A few weak-H/3 objects only 


TWFBG98 H/3 velocity dispersion correction 


+25% 


Weak-H/3 objects with high a only 


Metalhcity errors, [Z/H]: 






Theoretical response functions (TB95) 


0.10 dex 


Zeropoint 


Metal-hne fitting functions 


0.05 dex 


Zeropoint 


Unknown O/Mg ratio 


100%'' 


Scale 


Unknown C abundance 


0.05 dex 


Zeropoint 


Effect of [E/Fe] 7^0 on tracks 


0.05-0.1 dex 


Zeropoint 


Enhancement errors, [E/Fe]: 






Theoretical response functions 


20% 


Scale 


Metal-Une fitting functions 


0.10 dex 


Zeropoint 


Unknown C and O abundance 


0.01 dex 


Zeropoint 


Effect of [E/Fe] t^'O on tracks 


0.0 dex 





^Outer regions only. 

''Assumes that [O/H] is always solar regardless of [Z/H]. 



Note. — When signs are given, they are in the sense that the standard models are wrong by that 
amount. If no sign is given, the error could have either sign. 
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Fig. 10. — Two tests of the model dependences of our results, for the central r^/S aperture. In the top row, the results of Tantalo et al. (1998a; TCB98) are plotted 
as a function of our ages, metallicities, and enhancement ratios, inferred here using enhancement model 3 (C down, O up, labelled "O"), which has an elemental mix 
similar to that of Tantalo et al. We have shifted their results so that NGC 221 (M 32) has an age of 4 Gyr, [Z/H] = —0.1, and [E/Fe] = 0, as suggested by them. Note 
the disagreement between the inferred [Z/H] and [E/Fe] ([a/Fe]) distributions in the two studies. This disagreement is due mainly to their use of different response 
functions to correct the line indices for [E/Fe]^ 0; in particular, no correction to (Fe) was apphed. The bottom row shows the results of substituting the Bertelh et 
al. (1994; "Padova") isochrones for the Worthey (1994; "W94") isochrones in model 4; all other ingredients (fitting functions, response to non-solar-neighborhood 
abundance ratios) remain the iame. using enhancement model 4 ("CO"). Apart from slight zeropoint and slope differences, the two sets of isochrones give very 
similar results (see Appendix N). In all plots, the long-dashed Une is the line of equivalence for the two sets of results; the short-dashed line is a linear least-squares 
fit (iteratively rejecting 3-(t oiitliers). 



metallicities can be determined from (Fe) alone. In other 
words, the ages and metallicities of the G93 galaxies are effec- 
tively defined by the (Fe)-H/3 line-strength diagram alone in 
the TCB98 scheme, and the enhancement ratios [E/Fe] are de- 
fined by the offset of the galaxies in the M g2-H/3 line-strength 
diagram (scaled by some factor from the Borges et al. 1995 
fitting function for Mg2). By not correcting (Fe) upwards for 
Fe-deficiency, TCB98 underestimate the metallicites [Z/H] and 
overestimate the ages t and enhancement ratios [E/Fe]. This 
matches the behavior of residuals seen in Figure |T^. 

These systematic effects cause TCB98 to find a much nar- 
rower spread in in [Z/H] in the centers of the G93 galaxies 
than we do, and also a much wider spread in [E/Fe]. The 
narrow spread in [Z/H] prevents them from finding any age- 
metallicity relation, which is a major focus of our Paper II; 
conversely, the broad spread in [E/Fe] causes them to find a 
strong age-enhancement ratio relation, which we do not find in 
Paper II. Further discussion of these trends is reserved to future 
papers. However, it is clear that the adopted response functions 
can have far-reaching consequences for parameter correlation 
studies. 

6.2. Other authors 

Few other authors have fitted stellar population parameters 
to B almer line data. Kuntschner (1998; Kuntschner & Davies 
1998) has studied the line strengths of a complete, magnitude- 
limited {Mb < -17) sample of early-type galaxies in the For- 
nax cluster, split evenly between ellipticals and lenticulars. 
He derives stellar population parameters but does not correct 
his line strengths for non-solar abundance ratios. In general, 
Kuntschner finds old ages for ellipticals but a wide spread in 
the ages of SOs. His data are of excellent quality, and we lump 
them together with the G93 sample and analyze them in parallel 
in Paper II. 

Using moderate-5/A^ (~ 30) long-slit and fiber spectroscopy, 
J0rgensen (1999) has studied the line strengths of 115 early- 
type galaxies in Coma. Of these galaxies, 71 have measured 
Mg/? (Fe), and H/3(the last with typical errors of 0.22 A). How- 
ever, Trager (1997) has shown that errors in H/3 of this magni- 



tude (typical of the Lick/IDS galaxy sample; TWFBG98) se- 
riously compromise the determination of stellar population pa- 
rameters through correlated errors in age, metallicity, and en- 
hancement ratio; errors of ^ 0. 1 A in H/3 are required to deter- 
mine ages to 10% or better and to reduce the correlated errors to 
insignificant levels. Further consideration of Jorgensen's work 
is reserved to Paper II. 

Vazdekis et al. (1997) fit their own SSP-equivalent models to 
three early-type galaxies, including NGC 3379 and NGC 4472 
studied here. Their derived ages are about 50% larger than ours, 
for a variety of reasons. Although H/3is included in the suite of 
data fitted, it is only one among many features used. The resul- 
tant models significantly Mnc/er-predict their own H/3strengths, 
and matching them would yield ages as young or younger than 
we find. Their high ages (and low metallicities) seem to be 
driven by the very red near-IR colors of their models at high 
[Z/H], which in turn may stem from the cool giant-branch tips 
of the Padova isochrones used. Since giant-branch tempera- 
tures are still in flux, we prefer the Balmer lines, which are less 
sensitive to stellar evolution uncertainties. 

Fisher, Franx & Illingworth (1995) studied the Une strengths 
of nearby field ellipticals and brightest cluster galaxies (BCG). 
All of the nearby ellipticals (seven galaxies) were drawn from 
G93; the data are consistent, and we have therefore not added 
them to this series of papers. The BCG data are of slightly 
lower quality as the galaxies are more distant. Fisher et al. com- 
pare their line strengths to the W94 models — ignoring [E/Fe] 
variations — and generally find old (f ^ 10 Gyr) mean stellar 
populations in the centers. However, two of nine BCGs have 
H/3 strengths indicative of intermediate-age populations, NGC 
2329 (Abell 569) and NGC 7720 (Abell 2634). 

Jones & Worthey (1995) developed a novel H7 index that 
has lower sensitivity to metallicity, and therefore in principle 
better age discrimination. Using W94 models, they applied 
this index to the center of M 32 and determined an SSP age 
of f « 5-7 Gyr Our age for this object is only f = 3.0 ± 0.6 Gyr 
(with [Z/H] = 0.00 ±0.05 dex; formal errors only). Jones and 
Worthey also fitted other Balmer indices (including different 
versions of H7), which gave similarly low ages. They were not 
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able to identify a reason for the discrepancy. This disagreement 
among Balmer indices is an outstanding issue. 

Finally, we mention the results of Rose (1985, 1994) for a 
sample of 10 normal elliptical nuclei, 6 of which overlap with 
our sample, including M 32. Rose's spectra were taken around 
the 4000 A break, and he developed a large number of stel- 
lar population indicators in this wavelength region, including 
the Balmer index Ca II H+He/Ca II K, a sensitive measure of 
the presence of hot A and B stars, and Sr II/Fe I, a measure of 
the total dwarf-to-giant light. By balancing these and other in- 
dices. Rose found that there must be a substantial intermediate- 
temperature component of dwarf light in all of these galaxies, 
but that no more than 2% of 4000 A light could come from stars 
hotter than FO. He concluded that all 10 galaxies contained a 
significant component of intermediate-age main sequence stars. 
Our SSP ages range from 3 to 10 Gyr for the 6 galaxies in com- 
mon, consistent with these conclusions. 

7. SUMMARY 

We have presented central (r(,/8) and global {re/2) line 
strengths for the Gonzalez (1993) local elliptical galaxy sample. 
A method for deriving SSP-equivalent stellar population pa- 
rameters is presented using the models of Worthey (1994), sup- 
plemented by model-atmosphere line-strength response func- 
tions for non-solar element abundance ratios by Tripicco & Bell 
(1995). The resultant stellar population parameters broadly 
confirm the findings of G93 in showing a wide range of ages 
but a fairly narrow range of metallicities and enhancement ra- 
tios. Differences among galaxies in the sample are larger than 
radial differences within them. 

Four different models are considered with different patterns 
of element enhancement. The best-fitting model (model 4) 
has all elements enhanced or normal except for the Fe-peak 
(and Ca), which are depressed. The actual atomic abundance 
ratios of the so-called "enhanced" elements are in fact virtu- 
ally solar — it is really the Fe-peak elements that are depressed. 
Indeed, the TB95 response functions imply that the observed 
strengthening of M.gb is not due to an overabundance of Mg 
but to an M«t/erabundance of Fe (and Cr). It is shown that C 
must also belong to the enhanced group (i.e., it does not follow 
Fe, as sometimes assumed). Hence, a more accurate descrip- 
tion of elliptical galaxies is that they failed to make Fe-peak 
elements rather than that they made an overabundance of a- 



elements. The element enhancement pattern of ellipticals will 
be considered in more detail in a future paper 

Sources of error in the population parameters are consid- 
ered. Contamination of H/3 by hot stars such as horizontal 
branch stars and blue stragglers can cause small reductions in 
the measured ages of the oldest galaxies but cannot noticeably 
affect the strong H/3 lines, and thi is the dedu ced low ages, of 
young el lipticals (as also found by [Rose 1985i 1994 and G reg- 
gio 1997). Emission fill-in may increase the measured ages 
of a few, largely old galaxies, but the broad age distribution 
is unaffected by whether any emission corrections are made or 
not. Uncertainties in the theoretical tracks, index response func- 
tions, element enhancement patterns, and the Lick/IDS metal- 
licity scale all affect the absolute zero points of age, [Z/H], and 
[E/Fe] at the level of a few tens of percent or tenths of a dex — 
but not the relative age rankings among galaxies. 

Finally, we have compared our population parameters to 
those derived by TCB98, who apply a different modeling tech- 
nique to the G93 sample. Our values of [Z/H] and [E/Fe] cor- 
relate with theirs, but the slopes differ significantly from unity. 
This appears to stem from the use of different response func- 
tions; in particular, TCB98 do not correct (Fe) for the under- 
abundance of Fe. When this is allowed for, the two studies are 
consistent. 

Future papers will discuss the central stellar populations of 
the G93 sample in detail, correlations between stellar popula- 
tions and structural parameters, scaling relation of these local 
ellipticals in the context of stellar populations, and stellar pop- 
ulation gradients in elliptical galaxies. 



We thank Drs. M. Bolte, A. Bressan, D. Burstein, J. Dal- 
canton, G. lUingworth, D. Kelson, I. King, A. McWilliam, A. 
Renzini, M. Rich, M. Salaris, and A. Zabludoff for stimulating 
discussions. We thank especially the referee. Dr. J. Rose, for a 
careful and thorough reading of the manuscript which helped 
improve the final presentation. We are indebted to Drs. M. 
Tripicco and R. Bell for calculating response functions for the 
Lick/IDS indices, without which this work would not have been 
possible, and to Dr Tripicco for sending electronic versions 
of their tables. We also thank Dr. Salaris for sending us his 
and Dr Weiss's solar-metallicity, a-enhanced isochrones in ad- 
vance of publication. 



REFERENCES 



Arimoto, N. 1996, in From Stars to Galaxies: The Impact of Stellar Physics on 
Galaxy Evolution, ed. C. Leitherer, U. Fritze-von Alvensleben, & J. Huchra, 
A. S. R Conf. Ser., vol. 98 (San Francisco: ASP), p. 287 

Bailyn, C. D. 1995, ARA&A, 33, 133 

Baum, W. A. 1959, PASP, 71, 106 

BertelU, G., Bressan, A., Chiosi, C., Fagotto, F, & Nasi, E. 1994, A&AS, 106, 
275 

Bessel, M. S., Brett, J. M., Scholz, M., & Wood, P R. 1989, A&AS, 77, 1 
Bessel, M. S., Brett, J. M., Scholz, M., & Wood, P R. 1991, A&AS, 89, 335 
Borges, A. C, Idiart, T. R, de Freitas Pacheco, J. A., & Thevenin, F 1995, AJ, 
110, 2408 

Bressan, A., Chiosi, C, & Tantalo, R. 1996, A&A, 311, 425 
Burstein, D., Faber, S. M., Gaskell, C. M., & Krumm, N. 1984, ApJ, 287, 586 
Burstein, D., Bertola, F, Buson, L. M., Faber, S. M., & Lauer, T. R. 1988, ApJ, 
328, 440 

Carrasco, L., Buzzoni, A., Salsa, M., & Recillas-Cruz, E. 1996, in Fresh Views 

on Elliptical Galaxies, eds. A. Buzzoni, A. Renzini, & A. SeiTano, A. S. P. 

Conf. Ser., vol. 86 (San Francisco: ASP), p. 235 
Castro, S., Rich, R. M., McWilHam, A., Ho, L. C, Spinrad, H., Filippenko, A. 

v., & Bell, R. A. 1996, AJ, 111, 2439 
Cayrel de Strobel, G., Soubiran, C, Friel, E. D., Ralite, N., & Francois, P 1997, 

A&AS, 124, 299 

Chaboyer, B., Demarque, R, Keman, R J., & Krauss, L. M. 1998, ApJ, 494, 96 



Chariot, S. & Bruzual, G. 1991, ApJ, 367, 126 

Chariot, S., Worthey, G., & Bressan, A. 1996, ApJ, 457, 625 

Couch, W. J. & Sharpies, R. A. 1987, MNRAS, 229, 423 

Davies, R. L., Sadler, E. M., & Peletier, R. F 1993, MNRAS, 262, 650 

de Vaucouleurs, G., de Vaucouleurs, A., Corwin, H. G., Buta, R. J., Paturel, 

G., & Fouque, P. 1991, Third Reference Catalogue of Bright Galaxies (New 

York: Springer- Verlag) (RC3) 
Dressier, A. & Gunn, J. E. 1983, ApJ, 270, 7 

Edvardsson, B., Andersen, J., Gustafsson, B., Lambert, D. L., Nissen, P. E., & 

Tomkin, J. 1993, A&A, 275, 101 
Faber, S. M. 1972, A&A, 20, 361 
Faber, S. M. 1973, ApJ, 179, 731 
Faber, S. M. & Jackson, R. E. 1976, ApJ, 204, 668 

Faber, S. M., Friel, E. D., Burstein, D., & Gaskell, CM 1985, ApJS, 57, 711 
Faber, S. M., Wegner, G., Burstein, D., Davies, R. L., Dressier, A., Lynden-Bell, 

D., & Terlevich, R. J. 1989, ApJS, 69, 763 
Faber, S. M., Trager, S. C, Gonzalez, J. J., & Worthey, G. 1995, in lAU 

Symposium 164, Stellar Populations, eds. P. C. van der Kruit & G. Gilmore 

(Dordrecht: Kluwer), p. 249 
Ferraro, F R., Fusi Pecci, F, & Bellazzini, M. 1995, A&A, 294, 80 
Fisher, D., Franx, M., & lllingworth, G. 1995, ApJ, 448, 119 
Fisher, D., Franx, M., & lllingworth, G. 1996, ApJ, 459, 110 
Gonzalez, J. J. 1993, Ph.D. thesis. University of California, Santa Cruz (G93) 



30 



STELLAR POPULATIONS OF EARLY-TYPE GALAXIES. I. 



Gorgas, J., Faber, S. M., Burstein, D., Gonzalez, J. J., Courteau, S., & Prosser, 

C. 1993, ApJS, 86, 153 
Gottlieb, D. J. & Bell, R. A. 1971, A&A, 19, 434 
Goudfrooij, P. & Emsellem, E. 1996, A&A, 306, L45 

Gratton, R. G., Fusi Pecci, F, Can'etta, E., Clementini, G., Corsi, C. E., & 

Lattanzi, M. 1997, ApJ, 491, 749 
Green, E. M., Demarque, P., & King, C. R. 1987, The Revised Yale Isochrones 

and Luminosity Functions (New Haven: Yale University Observatory) 
Greggio, L. 1997, MNRAS, 285, 151 

Grevesse, N., Noels, A., & Sauval, A. J. 1996, in Cosmic Abundances, eds. S. 
S. Holt & G. Sonneborn, A.S.P Conf. Ser, vol. 99 (San Francisco: ASP), p. 
117 

Giillmair, C. J., et al. 1996, AJ, 112, 1975 

Grundahl, F, Vandenberg, D. A., & Andersen, M. I. 1998, ApJ, 500, L179 
Gustafsson, B., Kja;rgaard, P, & Anderson, S. 1974, A&A, 34, 99 
Gunn, J. E. & Stryker, L. L. 1983, ApJS, 52, 121 
Hansen, L. & Kja:rgaard, P 1971, A&A, 15, 123 

Ho, L. C, Filippenko, A. V., & Sai'gent, W. L. W. 1997, ApJS, 112, 315 

Iglesias, C. A.. Rogers, R J., & Wilson, B. G. 1992. ApJ, 397, 717 

Jensen, J. B., Luppino, G. A., & Tonry, J. L. 1996, ApJ, 468, 519 

Jensen, J. B., Tonry, J. L., & Luppino, G. A. 1998, ApJ, 505, 111 

Johnson, H. L. 1966, ARA&A, 4, 193 

Jones, L. A. & Worthey, G. 1995, ApJ, 446, L31 

J0rgensen, I. 1997, MNRAS, 288, 161 

J0rgensen, I. 1999, MNRAS, 306, 607 

Kuntschner, H. 1999, Ph.D. thesis. University of Durham 

Kuntschner, H. & Davies, R. L. 1998, MN, 295, L29. 

Kurucz, R. L. 1992, unpubhshed 

Lee, Y.-W. 1994, ApJ, 430, 113 

Liebert, J., Saffer, R. A., & Green, E. M. 1994, AJ, 107, 1408 

Matteucci, F 1994, A&A, 288, 57 

McClure, R. D. & van den Bergh, S. 1968, AJ, 73, 313 

McWilliam, A. 1997, ARA&A, 35, 503 

McWilHam, A. & Rich, R. M. 1994, ApJS, 91, 749 

O'Connell, R. W. 1976, ApJ, 206, 370 

O'Connell, R. W. 1980, ApJ, 236, 430 

Peletier, R. 1989. Ph.D. thesis, Rijksuniversiteit Groningen 

Pont, F, Mayor, M., Turon, C, & Vandenberg, D. A. 1998, A&A, 329, 87 

Rabin, D. 1982, ApJ, 261.85 

Reid, L N. 1997, AJ, 114, 161 

Reid, L N. 1998, AJ, 115,204 

Renzini., A. 1986, in Spectral Evolution of Galaxies, eds. C. Chiosi & A. 

Renzini (Dordrecht: Kluwer), p. 1 5 1 
Rich, R. M. 1988, AJ, 95, 828 



Rich, R. M., et al. 1997, ApJ, 484, L25 

Robinson, L. B. & Wampler, E. J. 1972, PASP, 84, 161 

Rose, J. A. 1985, AJ, 90, 1927 

Rose, J. A. 1994, AJ, 107, 206 

Salai-is, M. , Chieffi, A. & Straniero, O. 1993, ApJ, 414, 580 
Salai-is, M. & Weiss, A. 1998, A&A, 335, 943 

Sandage. A. & Bedke, J. 1994, The Carnegie Atlas of Galaxies (Washington, 

D.C.: Carnegie Institution of Washington with the Flintridge Foundation) 
Sandage. A. & Tammann, G. A. 1987, A Revised Shapley-Ames Catalog 

of Bright Galaxies, 2nd ed. (Washington, D. C: Carnegie Institution of 

Washington) (RSA) 
Schweizer, S., Seitzer, P, Faber, S. M., Burstein, D., Dalle Ore, C. M., & 

Gonzalez, J. J. 1990, ApJ, 364, L33 
Schweizer, F. & Seitzer, P 1992, AJ, 104, 1039 
Spinrad, H. & Taylor, B. J. 1971, ApJS, 22, 445 

Sosin, C, et al. 1997, in Advances in Stellai' Evolution, eds. R. T. Rood & A. 

Renzini (Cambridge: Cambridge University Press), p. 92 
Tantalo, R., Chiosi, C, & Bressan, A. 1998, A&A, 333, 419 (TCB98) 
Tantalo, R., Chiosi, C, Bressan, A., Marigo, R, & Portinari, L. 1998b, A&A, 

335, 823 

Trager, S. C. 1997, Ph.D. thesis. University of Cahfornia, Santa Cruz 

Trager, S. C, Worthey, G., Faber, S. M., Burstein. D., & Gonzalez, J. J. 1998, 

ApJS, 116, 1 (TWFBG98) 
Tripicco, M. & Bell, R. A. 1995, AJ, 110, 3035 (TB95) 
VandenBerg, D. A. 1985, ApJS, 58, 711 
VandenBerg, D. A. & Bell, R. A. 1985, ApJS, 58, 561 
VandenBerg, D. A. & Laskarides, P G. 1987, ApJS, 64, 103 
Vazdekis, A., Casuso, E., Peletier, R. F, & Beckman, J. E. 1996, ApJS, 106, 

307 

Vazdekis, A., Peletier, R. F, Beckman, J. E., & Casuso, E. 1997, ApJS, 111, 
203 

Weiss, A., Peletier, R. R, & Matteucci, R 1995, A&A, 296, 73 

Wheeler, J. C, Sneden, C. & Tmran, J. W. 1989, ARA&A, 27, 279 

Woosley, S. E. & Weaver, T. A. 1995, ApJS, 101, 181 

Worthey, G. 1992, Ph.D. thesis. University of California. Santa Cruz 

Worthey, G. 1993, ApJ, 415, L91 

Worthey, G. 1994, ApJS, 95, 107 (W94) 

Worthey, G. 1998, PASP, 110, 888 

Worthey, G., Faber, S. M., & Gonzalez, J. J. 1992, ApJ, 398, 69 

Worthey, G., Faber, S. M., Gonzalez, J. J., & Burstein, D. 1994, ApJS, 94, 687 

Worthey, G. & Ottaviani, D. L. 1997, ApJS, 111, 377 

Yi, S., Lee, Y.-W., Woo, J.-H., Park, J.-H.. Demarque. P & Oemler, A., Jr 1999, 
ApJ, 513, 128 



APPENDIX 

STELLAR POPULATION PARAMETERS USING PADOVA ISOCHRONES 

This section provides more details on the "Padova" models discussed in Section |6.1| . The Padova models are identical to the 
W94 models except that the isochrones (and opacities) are replaced by the isochrone library of Bertelli et al. (1994). This isochrone 
library is based on the stellar evolutionary tracks developed by the Padova group (see Bertelli et al. 1994 and Chariot, Worthey, & 
Bressan 1996 for more details) using the Iglesias, Rogers & Wilson (1992) radiative opacities. These isochrones include all phases 
of stellar evolution from the ZAMS to the remnant stage for stars of masses in the age range 0.004 < f < 16 Gyr and metallicity 
range 0.0004 < Z < 0.1 (Zq = 0.02). The models include convective overshooting in stars more massive than 1 Mq and an analytic 
prescription for the TP-AGB regime. 



Figure Al presents the inferred H/3, Mgfe, and (Fe) line strengths for the W94 models using the Padova isochrones. Comparing 
this with Figure |l] shows small differences: the Padova models have a higher metallicity at a given Mgb or (Fe) strength and a 
younger age at a given H/3 strength. 

Figure HQ (bottom row) shows the results of applying these new models to the G93 central (r(,/8) line strengths. The derived ages, 
metallicities, and enhancement ratios agree quite well with the W94 models, apart from slight slope changes and zeropoint offsets: 

lOgfpadova = 1-02 logfw94-0.10, (Al) 

[Z/H]padova = 1-09 [Z/H]^,4 + 0.03, (A2) 
[E/Fe]p3,_ = 0.97 [E/Fe]^,^ + 0.02. (A3) 

The above are linear least-square fits using enhancement model 4 and rejecting 3-cr outliers. The fit for log? is in accordance with 
the results of Chariot, Worthey & Bressan (1996): changing isochrones can alter the inferred ages from line strengths at young ages 
by as much as ^ 25%; agreement at old ages is within 10%. On average, the inferred metallicities [Z/H] are increased by « 10% in 
the Padova models, as expected at fixed line strengths from the 3/2 rule (Alogf/A[Z/H] « 1.4 between the two sets of models). 

CORRECTIONS TO H/3 



Emission corrections 
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Padova models 
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Fig. A1. — The effect of substituting model isoclirones from Bertelli et al. (1994) for those of Worthey (1994). The grid in this figure should be compared to that 
in Figure 1. The basic morphology of the grid is unchanged. The data are the central r^/S values as in Figure 1. 
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T his se ction presents stellar population parameters derived by omitting the emission fill-in correction to H/3 discussed in Sec- 
tion 2.2.2. This follows the suggestion by Carrasco et al. (1996) that no correction to li/3 should be made for residual H/? emission 
based on fO III], as they find no such correlation in their own sample of early-type galaxies. 

Figure Bl presents the Balmer-metal-line diagrams for the G93 galaxies through the r^/Bi aperture with the H/3 correction omitted. 
As expected from the additive nature of the correction, galaxies now appear lower in the grid, and therefore older and more metal- 
poor than in Figure |l[ Neglecting emission corrections forces some galaxies to have unreasonably old ages: for example, without 
corrections, NGC 1453, NGC 2778, NGC 4261, NGC 4374, NGC 5813, NGC 5846, NGC 7052 have ages ^ 20 Gyr. Since all 
of these galaxies have clear H/3 emission (see Figs. 3.11 and 4.10 of G93), omitting the corrections makes no sense. Furthermore, 
careful checking reveals that a few galaxies (e.g., NGC 4552, NGC 4649, NGC 5813, NGC 5846, NGC 7052) actually appear to 
have H/3 a little stronger than the standard ratio and are therefore probably undercorrected in our standard treatment. Fixing them 
would move them up by a few hundredths in log H/3 and decrease their ages by <^2Q%. Since some of these are also objects that lie 
low in the grid, this correction would improve their positions. 




Fig. Tj ' — - Line strengths of early-type galaxies in tlie Gpnzalez (1993) sample in the central r^/S aperture, but now with H/3 emission corrections omitted 
(Section |2.2.2| ). W94 model grids are overlaid as in Figure hi Galaxies on average lie a little lower here than in Figure 1 but their relative parameters (including 
ages) are little affected, showing that SSP pai'ameters are rather insensitive to the exact H/3 correction used. 
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Figure ^T] without H/? corrections looks essentially like the original one — the large age spread and relative parameter rankings of 
the galaxies are essentially the same. This point is reinforced in the histograms of Figure B2, which are nearly identical to those in 
Figure ^ We conclude that emission corrections are needed to derive the best age estimates for early-type galaxies, but that their 
exact magnitude does not affect our broad conclusions. 

Velocity dispersion corrections 

This section derives a third set of population parameters using the Lick/IDS velocity dispersion corrections of TWFBG98 for H/3 
rather than the template-based corrections of G93. As G93 does not provide raw H/3 strengths, we use his Figure 4.1 to estimate his 
velocity corrections and use the m to " uncorrect" the fully corrected line strengths back to raw strengths (after removing the emission 



correction discussed in Section 2.2.2). For most galaxies, G93's velocity corrections are insignificant, but for high-cr galaxies they 



tend to be negative. After the G93 corrections are removed, we a pply t he positive corrections presented in TWFBG98 (their Figure 3) 



and then reapply the emission corrections discussed in Section [2.2.2. The newly corrected line strengths are plotted in Figure B3 



and the resulting stellar population parameter histograms are presented in Figure B4 

The new corrections move only a few high-cr galaxies, and therefore the relative age ranking of the sample is unaffected. The 
affected galaxies again tend to lie at the bottom of the grid and are again moved up by the new corrections (by ^0.04 in log H/3) so 
that galaxies that formerly lay below the grid at high ages now tend to lie on it. This correction, like the refined emission corrections 
of the previous section, thus improves the ages of the oldest objects. 

The G93 H/3 velocity corrections were based on a very high-S/N stellar template fit to each galaxy, whereas the TWFBG98 
corrections are based on a statistical average over stellar spectral types whose correction curves scattered widely. Nevertheless, it 
is possible that the TWFBG98 corrections are actually more acurate. The G93 stellar templates are a superb match to most of the 
spectrum except at H/3, where emission corrupted the data. Hence the template match (and correction) at H/3 in particular may be 
poor. The TWFBG98 corrections were selected to match a large number of stars with about the same H/3 strength as typical galaxies 
and could therefore be better on average. 
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Fig. B2. — The effect of H/3 emission corrections on the SSP parameter distiihiitions of the Conzalez (1993) sample using enhancement model 4 (central /8 
aperture). Solid histograms use the H/3 emission fill-in correction (Section [2.2.2| ; see Figure uk: dotted histograms do not include this correction. Without the 
emission correction, galaxies are typically slightly more metal-poor and older (NGC 1453 is '^26 Gyr old in this model), but their relative ranking in the various 
parameters is little affected. 
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Lick/IDS H(S velocity dispersion corrections 
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Fig. B3. — Line strengths_aL£arly-type galaxies in the Gonzalez (1993) sample in the central rc/S aperture but naa/ using the Lick/IDS velocity dispersion 
corrections for H/3 (Section p.2.l| ) rather than those of G93. The oldest galaxies lie sUghtly higher here than in Figures hi and H and their SSP-equivalent ages are 
reduced by about 25% (see Figure 16). Overly old ages of low-lying galaxies through the re/2 aperture are similarly reduced (cf. Figure 0). 
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Lick/IDS velocity dispersion corrections 
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Fig. B4. — The effect of H/3 velocity corrections on the SSP parameter distributions nf the Gonzalez (1993) sample using enhancement model 4 (central f>/8 
aperture). Solid histograms use the G93 H/3 velocity dispersion correction (Section |2.2.l| ; see Figure bl); dotted histograms use the Lick/IDS velocity dispersion 
correction. With the Lick/IDS correction, galaxies are younger, slightly more metal-rich, and slightly more enhanced in a-elements, but retain their basic relative 
parameter rankings. 



